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Large-signal  analysis  programs  are  widely  used  by  the  microwave 
tube  community  to  predict  the  performance  of  helix  traveling-wave 
tubes.  These  nonlinear  programs  allow  an  engineer  to  optimize  a  tube 
design  without  building  several  prototype  tubes.  Unfortunately,  cur¬ 
rently  available  programs  are  complex  and  poorly  documented.  The  pur¬ 
pose  of  this  project  Is  to  provide  a  structured,  well  documented  large- 
signal  analysis  program  to  the  scientific  and  industrial  community. 

The  program,  TWT1D,  is  based  on  the  one-diraenslon  nonlinear  theory 
developed  by  M.  K.  Scherba  and  J.  E.  Rowe.1  The  program  was  written  in 
FORTRAN77  and  implemented  on  an  HP-1000  computer  at  Teledyne  MEC  in  Palo 
Alto,  California.  Complete  documentation  for  the  program.  Including  a 
source  code  listing  and  a  User's  Manual,  is  given  in  this  report. 
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II.  LARGE-SIGNAL  THEORY 


The  large-signal  theory  of  traveling-wave  amplifiers  differs  from 
small-signal  theory  in  several  ways.  Perhaps  the  most  significant 
difference  is  the  way  each  theory  models  the  electron  beam.  In  small- 
signal  theory,  the  beam  is  modeled  as  drifting  charged  fluid.  Its 
charge  density  and  electron  velocity  are  single  valued  functions  of 
distance,  and  electrons  cannot  overtake  one  another.  This  treatment  of 
the  electron  beam  is  called  an  Eulerian  formulation.  By  comparison,  the 
Lagrangian  analysis  used  in  large-signal  theory  breaks  the  beam  up  into 
representative  charge  groups.  This  approach  allows  charge  density  and 
electron  velocity  to  be  multivalued  functions,  and  electrons  to  overtake 
one  another.  Another  difference  between  the  theories  is  the  way  they 
treat  the  force  and  continuity  equations.  The  equations  are  linearized 
in  small-signal  theory,  so  the  theory  is  limited  to  circuits  with  small 
"c"  and  small  input  signals.  Linearizing  the  equations  allows  them  to 
be  solved  analytically.  The  solution  is  a  single  number  which  gives  the 
gain/unit  length  for  the  tube.  By  contrast,  the  large-signal  theories 
have  no  analytic  solution  and  must  be  solved  on  a  computer.  The  result¬ 
ing  solution  is  a  profile  of  power  versus  axial  position. 

The  large-signal  theory  offers  many  advantages  over  small-signal 
theory.  Since  the  equations  are  nonlinear,  harmonics  can  be  generated 
and  followed  down  the  tube.  This  allows  the  effects  of  intermodulation 
to  be  studied.  Unlike  small-signal  theory,  large-signal  theory  can  be 
used  to  predict  the  saturation  power  and  efficiency  of  a  tube.  In 
addition,  using  a  Lagrangian  analysis  means  that  electron  trajectories 
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can  be  followed  down  the  tube  so  beam  bunching  can  be  observed.  Final¬ 


ly,  the  theory  can  accommodate  attenuators,  severs,  and  velocity 


changes. 


The  first  attempt  at  developing  a  large-signal  theory  for  helix 


traveling-wave  tubes  was  made  by  A.  T.  Nordsieck2  in  1953.  He  was  the 


first  to  use  Lagrangian  analysis  to  model  the  interaction  between  the 


circuit  and  the  electron  beam  for  a  single  input  signal.  However,  his 


theory  was  subject  to  several  limitations.  To  begin  with,  he  assumed 


that  space-charge  forces  could  be  neglected.  Secondly,  he  limited  his 


theory  to  lossless  circuits  with  small  C.  A  year  later,  H.  C.  Poulter3 


extended  the  theory  to  include  space-charge  forces  and  to  account  for 


circuit  loss  and  finite  C.  Two  years  after  that,  J.  E.  Rowe1*  amended 


the  theory  to  allow  for  circuits  with  a  large  C.  Finally,  in  1971,  M. 


K.  Scherba  and  J.  E.  Rowe1  expanded  the  theory  to  include  multiple  input 


signals.  It  is  this  theory  that  is  the  basis  for  the  program  TWT1D. 


In  Rowe's  model,  the  helix  is  represented  by  a  lumped  element 


transmission  line,  and  the  electron  beam  is  represented  by  disks  of 


charge.  These  charge  groups  move  along  the  circuit  inducing  charges  in 


the  circuit.  The  RF  field  in  the  circuit  interacts  with  the  electron 


disks  and  causes  them  to  bunch.  This  model  was  first  used  by 


Brillouin,5  and  is  schematically  represented  in  Fig.  1. 


IS 


m 


Fig.  1.  Schematic  model  of  a  helix  with  an  electron  beam 
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inhere 


dg  =*  length  of  section 

Rq  -  Rdg  »  resistance/section 

Lq  *  LdQ  *  inductance/section 

Cq  *  CdQ  *  capacitance/section 

In  **  current  in  Lq  of  section  n 

Vn  =  potential  on  capacitor  n 

Qn  *■  CqVo  ■  charge  on  capacitor  n 

q*  »>  p  d_»  charge  in  the  beam  in  section  n 
n  n  0 

The  circuit  equation  for  this  transmission  line  is 


£i.R9V  2  3i  (ift  R  3 A 

U2  +  LH  »0}t2  V0^2 


where  the  following  definitions  have  been  used: 


Vq  -  /l/LC  =■  the  characteristic  phase  velocity  of  the  circuit 


To  express  this 
used : 


■  the  characteristic  impedance  of  the  circuit 
equation  in  Pierce6  notation,  the  following  relation  is 


R/L  =  2usCd 


(2) 


where  C  is  Pierce’s  gain  parameter  and  d  is  his  loss  parameter.  The 
equation  then  becomes 


32V 

3z2 


+  2u>Cd 


sv  „  7  ( lie. 

3t "  °  °Lt2 


2u)Cd  It 

3  t  i 


(3) 


For  multiple  signals,  there  will  be  a  circuit  equation  for  each  signal: 


3V_ 

+  2u  C  d  x-z- 
n  n  n  dt 


n 


v0nZ0n' 


o  +  2o)  C  d 
2  n  n  n 


(4) 


A.  Phase  Relation  Equation 

In  his  Lagrangian  analysis,  Nordsieck  introduced  two  independent 
variables.  They  are  the  normalized  distance,  y,  and  the  entry  phase  of 
the  fundamental  component  of  the  input  signal,  (fr^.  They  are  defined  as 


C  u> 

A  1  1 

y  -  -  Z 

uo 


U)  Z 

•  ft-i - «  t 

°i  uo  10 


(5) 


(6) 
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where  is  the  fundamental  frequency  in  the  radians,  and  Cj  is  the  gain 
parameter  for  the  fundamental  frequency.  uq  is  the  dc  stream  velocity 
in  m/a,  and  t  is  the  axial  position  on  the  tube  measured  in  meters. 
Using  these  new  variables,  the  circuit  equation  becomes 


.2  2  fvOn^  .  .  _  .  3Vn(y,t) 

Cl"l  (iH  ~2 +  Wn ~t 

\  0  /  3y 


a  p  (y,t)  ap  (y,t)' 

v0n20n  I  ~ -  7  2  *  2W*n  "TT 

at 


(7) 


Nordsieck  also  introduced  the  following  dependent  variables: 


♦>•‘01)  *  ■ c)  ■  9»<?)  ■  ^ y  ■  v  ■  vy) 


(8) 


a 

Z  dt 


-  C7&7  H1  -  “0[1  * 


(9) 


t,z  ri 


where  *  is  the  instantaneous  charge  group  phase  in  radians,  and  6  is 
n  n 

the  phase  lag  in  radians  of  the  circuit  wave  relative  to  a  traveling 
wave  with  phase  velocity  uQ,  This  phase  lag  is  due  to  beam  loading  as 
the  wave  gets  energy  from  the  beam.  uz  is  the  electron  velocity  normal¬ 
ized  to  the  initial  average  electron  velocity,  Uq.  Taking  the  deriva¬ 
tive  with  respect  to  time  of  Eq.  8  and  substituting  Eq.  9  yields 


aen(y>  a*n(r»*0l)  _  J^n_ 

5y  5y  "  u1C1 


1  *  2Cl»l*-«0l 


T 


(10) 


This  is  the  first  of  Rowe's  large-signal  equations.  It  is  the  velocity- 
phase  equation  relating  uz,  0n»  and 


B.  Circuit  Wave  Equations 


The  next  two  large-signal  equations  arise  from  the  circuit  equa¬ 


tion.  To  evaluate  the  circuit  equation,  define  the  RF  voltage. 


where 


v  (y.*01)  "  E  v(y>  v(*n) 

n  u 


ZmIn  -j*_ 

— — -  A  (y)  and  V(#  )  *  e 

n  n 


Z0i  is  the  interaction  impedance  for  the  fundamental  frequency,  and 
A  (y)  is  the  normalized  circuit  voltage  amplitude  for  the  nth  fre¬ 
quency.  Substituting  the  derivataives  of  Eq.  11  into  Eq.  4  yields  an 
equation  in  terms  of  sin  and  cos  Since  these  terms  are  mutually 
orthogonal,  they  can  be  separated  and  the  original  circuit  equation  can 
be  expressed  as  two  equations. 


d^A  (y)  / u>  d9  \  ^ 


+  Cnbn\2 


Rr  cos  (p)  /uq 


“l^O^O  '  0n 


An(y) 


2  2/  2 

d0  /»  \  /l  +  C  b  \ 

n  <1/,  j  i  n  W  n  n  \ 

dy2  ^){  c,  ) 


<Vr)  (  “n  d9„ 


dy  V1C1  "  dy_ 


R  sin  (p)  (u 
n  /  0 


❖  iVo  VO" 


where  R„  sin  (p  )  and  R  cos  (p)  are  the  sine  and  cosine  parts  of  the 
n  n 

space-charge  expression,  and  bn  is  Pierce's  velocity  parameter.  Now  the 
problem  remains  of  expressing  the  space-charge  expression  in  terms  of 
sin  #  and  cos  ♦  .  Since  the  beam-charge  density  is  rich  in  harmonics, 
it  is  convenient  to  expand  it  into  a  Fourier  series: 


*  sin  m*  2*  «  cos  2* 

p  (y»*  )  =»  J  - 2-  /  p  sin  ra$  d$  +5"  - —  f  p  cos  m*  d* 

rr  n'  ,  if  n  n  L  if  J  _  n 

ra=l  0  m=l  0 


The  p  inside  the  integral  may  be  determined  from  continuity  argu¬ 
ments  : 


1  +  2CiUly’\)l 


Using  this  expression  inside  the  Fourier  integral  yields 

"  '  “o’  l  “  "  0  1  +  2C 


2f  COS  m*  (y,*'  )d»' 

.  r  n  0101 

+  cos  m*  J  - .  —=7! j - t-j — i— 

n  1  +  ZCjuU.^iJ  _ 


The  summation  over  m  represents  the  sum  of  all  harmonic  frequencies  for 
each  of  the  n  fundamental  frequencies.  Assuming  that  the  circuit  beam 
interactions  at  the  harmonic  frequencies  are  negligible  and  keeping  only 
the  m  -  1  terms, 
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C.  Electron  Force  Equation 

Rove'*  final  large-signal  equation  is  derived  from  the  Lorentz 
force  equation. 


(21) 


where  and  V  are  the  space-charge  and  circuit  terns,  respectively. 

8  C  C 

Using  Eq.  9  and  the  chain  rule  for  derivatives  and  substituting  into  Eq. 
21  yields 
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where  Egc  la  the  space-charge  field.  The  constant  in  front  of  the 
circuit  wave  term  can  be  simplified  using  the  definition  of  Pierce's 
gain  parameter. 


a“(y,#n, ) 

-  [l  ♦  2ClUCy,#01 )]  -  C,  1 
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Rowe**  derived  a  space-charge  expression. 


where  F(#  -  )  is  a  space-charge  weighting  function  derived  by 

Foul ter. 5  Rowe  developed  a  cloeed  fora  expression  for  this  function: 


-  (•  -  •A,) 


.In  (•  -  ) 


♦  ten 


nFbTfuT7b1 


(»  -  *ni) 


for  0  <  (*  -  «’  )  <  2r 


where  s'  end  b'  ere  the  helix  and  beaa  radii,  respectively,  and 
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at  constant  a'/b' 


R„  Is  the  electron  plasma  frequency  reduction  factor, 
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where  ln  and  are  modified  Bessel  functions.  Substituting  the  space- 
charge  expression  Into  Eq.  22  yields  the  fourth  and  final  large-signal 
equation. 
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D.  Initial  Conditions 


The  solution  ot  Rowe's  large-signal  equations  may  be  treated  as  an 
i n i 1 1  a  1 -va  1  ue  problem  It  It  Is  assumed  that  the  helix  is  terminated  In 
its  characteristic  Impedance  so  that  there  .ire  no  reflections  at  the 
output.  The  initial  conditions  for  the  beam  and  circuit  wave  are 
derived  in  this  section.  In  deriving  these  conditions,  it  was  assumed 
that  the  electron  beam  is  unmodulated  at  the  input. 

At  the  input  to  the  tube  (v  =  *>),  there  is  no  RF  interaction,  so 
the  second  derivatives  are  zero.  Also,  since  the  beam  is  unmodulated  at 
the  input,  the  space-charge  terms  are  zero.  Thus,  the  circuit  wave 
phase  equation  becomes 


d  6  w  C  b 
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Substituting  this  into  the  circuit  wave  amplitude  equation  at  the  input 
vields 
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(  30 ) 


The  unmodulated  beam  at  the  input  is  represented  hv  uniformly 
distributing  the  electron  iharge  groups  over  one  cv<  le  of  the  fundamen¬ 
tal  frequency.  For  1  iharge  groups,  the  initial  phases  are  given  by 


Virtually  all  helix  traveling-wave  tubes  use  attenuators  and/or 
severs  to  prevent  oscillations.  Therefore,  any  program  used  for  the 


design  of  these  tubes  must  be  able  to  model  circuit  loss. 

The  simplest  way  to  model  loss  in  a  traveling-wave  tube  is  to 
increase  Pierce's  loss  parameter,  d.  This  is  the  method  recommended  by 
Rowe.7  Figure  2  shows  the  result  of  various  Gaussian  shaped  attenuators 
on  output  power,  and  Fig.  3  shows  the  result  for  a  region  of  uniform 
loss.  As  we  would  expect,  the  attenuators  increase  the  saturation 
length  without  decreasing  the  saturation  power.  (This  would  not  be  the 
case  if  the  attenuator  was  too  close  to  the  output.) 

Unfortunately,  there  is  a  problem  with  this  method  of  modeling 
loss  in  a  circuit.  It  will  lead  to  instabilities  for  large  values  of 
loss.  This  is  due  to  the  fact  that  both  the  propagation  constant,  6, 
and  the  impedance,  Z,  vary  with  the  resistance  in  the  circuit.  The 
change  in  6  causes  the  circuit  wave  to  lose  synchronism  with  the  beam, 
and  the  change  in  Z  causes  reflections  in  the  attenuator  due  to  the 
impedance  mismatch.  The  result  of  these  problems  is  that  the  circuit 
wave  phase  angle  oscillates  in  a  region  of  heavy  attenuation.  This  is 
illustrated  in  Fig.  4.  The  higher  the  loss  and/or  the  longer  the  region 
of  attenuation,  the  greater  the  oscillations.  In  order  to  work  around 
this  problem  in  the  model,  attenuators  are  replaced  by  severs  when  their 
1  os 8  exceeds  a  certain  level. 
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In  fact,  the  effect  of  a  strong  attenuator  Is  much  the  same  as 
that  of  a  sever.  In  a  sever,  the  circuit  wave  is  eliminated  and  the 
beam  drifts  down  the  tube  affected  only  by  the  space-charge  forces. 
Similarly,  a  strong  attenuator  will  diminish  the  circuit  wave  so  much 
that  the  electron  beam  is  unaffected  by  it.  Thus,  it  is  reasonable  to 
model  a  region  of  strong  attenuation  with  a  sever. 

The  equations  inside  the  sever  are  simply  Rowe's  large-signal 
equations  with  the  circuit  wave  terms  removed, 
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The  first  equation  is  the  electron  force  equation,  and  the  second  is  the 
electron  phase  equation.  Once  inside  a  sever,  the  circuit  wave  phase, 
the  circuit  wave  amplitude,  and  their  derivatives  are  zero.  The  remain¬ 
ing  two  drift  tube  equations  are  integrated  through  the  sever. 

At  the  end  of  the  sever,  the  circuit  wave  must  be  restarted. 
Rowe7  suggests  that  the  following  initial  conditions  be  used:  A^y)  * 

0.  9  (y)  *  0,  and  d0  (y)/dy  ■  -  w  b  c  /u.C,.  However,  the  last  condi- 
n  n  nnnll 

tion  requires  the  beam  to  be  unmodulated  (Eq.  28).  In  fact,  the  beam  is 
highly  modulated  in  the  sever,  and  it  is  this  modulation  that  restarts 
the  circuit  wave.  I  suggest  the  following  conditions  for  the  end  of  a 
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sever.  Since  there  Is  no  RF  interaction  in  the  sever,  the  second  deriv¬ 
atives  are  zero.  So  the  circuit  wave  phase  equation  becomes 
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and  the  circuit  wave  amplitude  equation  becomes 
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Note  that  Eq.  35  has  the  circuit  wave  amplitude,  An(y),  in  the 
denominator.  This  prevents  us  from  assigning  a  value  of  zero  to  this 
variable  at  the  end  of  a  sever.  Instead  of  zero,  the  maximum  of  1 
percent  of  the  input  value  or  1  percent  of  the  present  value  of  A^  is 
used.  This  value  is  large  enough  to  prevent  numerical  overflow,  yet 
small  enough  to  realistically  model  a  sever.  The  power  and  circuit  wave 
phase  are  plotted  using  both  Rowe's  and  the  new  initial  conditions  at 
the  end  of  a  sever  in  Figs.  5  and  6.  As  the  figures  show,  the  new 
initial  conditions  greatly  reduce  the  oscillations  in  the  circuit  phase 
plot  and  eliminate  them  completely  in  the  power  plot. 
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IV.  COMPUTER  IMPLEMENTATION  OF  ROWE'S  EQUATIONS 


This  section  describes  how  the  large-signal  equations  derived  in 
the  previous  sections  are  set  up  and  solved  in  TWT1D.  In  addition,  it 
shows  how  the  program  is  organized,  and  how  attenuators,  severs,  veloc¬ 
ity  steps,  and  tapers  are  handled.  TWT1D  is  composed  of  a  main  program 
block  and  23  subprograms.  The  program  blocks  can  be  grouped  according 
to  the  functions  they  serve,  as  illustrated  in  Table  1. 

Table  1.  TWT1D  program  blocks  organized  by  function. 


Function 

Program  Block 

PROGRAM  CONTROL 

TWT1D 

ZSORT 

EVENT 

INTSUB 

ROWE'S  LARGE  SIGNAL 

DIFFEQ 

EQUATIONS 

DIFF 

DRFTEQ 

DRIFT 

INTEGRALS  IN  ROWE'S 

FINT 

EQUATIONS 

INTI 

INT2 

NUMERICAL  ANALYSIS 

DESOLV 

STEP 

INTRP 

MACHIN 

BESSEL 

OUTPUT 

OUTPUT 

OUTRS 

CIRCUIT  CHANGES 

ATTEN 

CHANGE 

DIAGNOSTICS 

ENERGY 

BADIN 

ROWE'S  SPACE-CHARGE 

SPACE 

EXPRESSION 

SIMINT 

A.  Program  Control 

TWTlD  is  set  up  to  Integrate  Rowe's  large-signal  equations  inwn  a 


tube,  with  stops  for  attenuators,  velocity  changes,  and  restart  data 
printouts.  At  each  of  these  stops,  the  appropriate  subroutine  is  ci’led 
to  implement  the  change  in  parameters  or  the  printout.  The  flow  of  the 
program  is  controlled  by  the  main  program  block,  and  subroutines  /.SORT, 
EVENT,  and  INTSU8. 

The  main  program  block  is  divided  into  9  sections.  The  first 
section  sets  up  the  arrays,  variables,  and  constants  used  in  the  main 
program  and  various  subprograms.  The  second  section  controls  the  ter¬ 
minal  display  and  keyboard  input  when  the  program  is  run.  The  third 
section  reads  the  data  from  the  input  data  file.  In  addition,  the 
attenuators  and  velocity  tapers  are  set  up  in  this  section.  The  next 
section  converts  the  Input  data  Into  Pierce  parameters  and  normalized 
Rowe  variables.  Section  5  reads  the  restart  data  or,  if  no  restart  data 
are  available,  Initializes  the  beam  and  circuit  wave.  In  Section  f> ,  the 
output  data  file  is  opened,  and  the  Pierce  parameters  computed  in  Sec¬ 
tion  4  are  printed  out.  In  order  to  facilitate  the  program  flow.  Sec¬ 
tion  7  puts  all  events  that  require  stopping  the  integration  into  a 
single  list,  sorted  according  to  their  position  along  the  circuit. 
Thene  events  are  restart  data  printouts,  attenuators,  and  velocity  steps 
and  tapers.  The  eighth  section  controls  the  actual  flow  of  the  pro¬ 
gram.  It  moves  down  the  axial  position  list,  integrating  Rowe's  large- 
signal  equations.  At  each  event,  the  integration  is  stopped,  the  type 
of  event  is  determined,  and  the  appropriate  subroutine  called.  The  last 


sect  ton  c  imputes  and  prints  out  the  saturation  power,  gain,  efficiency, 
iii  length.  The  output  files  are  then  closed,  and  program  execution 
st  ips . 

Subroutine  TS'iRT  does  the  actual  sorting  of  the  axial  position 
list.  It  creates  two  Lists:  one  of  stop  positions,  ZSTOP(l),  and  one 
it  event  tvpe  c  i  ies  corresponding  to  the  stop  positions,  SWLIST(l). 
Sach  event  tvpe  Is  assigned  a  specific  prime  number  code.  When  more 
than  me  event  occurs  it  a  position,  the  event  type  codes  are  multi¬ 
plied.  Thus,  the  lumbers  can  he  factored  to  determine  which  events 
occur  it  each  st  ip  position. 

Subroutine  F1VKNT  checks  the  event  code  at  each  stop  position  to 
deter  nine  which  events  occur.  The  event  codes  used  are:  2  Print 
restart  lata;  1  Start  attenuator;  S  Stop  attenuator;  7  velocity  change. 

The  control  code  In  the  main  program  block  stops  the  Integration 
it  the  large-signal  equations  for  attenuators,  velocity  changes,  and 
restart  data  printouts.  however,  the  Integration  must  also  be  stopped 
fir  the  output  printouts.  Subroutine  INTSUB  accomplishes  this  by  break¬ 
ing  up  tube  sections  between  successive  stop  events  into  print  Inter¬ 
vals.  \t  each  print  Interval,  subroutines  are  called  to  Integrate 
Rowe's  equations  and  print  the  output  data.  In  addition,  INTSUB  ensures 
that  the  equations  are  Integrated  to  the  end  of  the  section. 

B.  Rowe's  Large-Signal  Equations 

Rowe's  circuit  wave  amplitude  and  circuit  wave  phase  equations  are 


second  order  differential  equations.  However,  the  differential  equation 


solver  used  la  TVT1D  solves  on  1  y  first  order  differential  equations. 
Moreover,  the  equations  must  be  of  the  form, 
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So,  the  equations  must  be  manipulated  to  fill  these  requirements.  Let 
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Substituting  these  new  variables  into  Rowe's  equations  and  rearranging 
yields 
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These  are  the  differential  equations  that  are  set  up  aad  evaluated  la 
subroutines  DIPPF.Q,  DIFF,  DRFTFQ,  and  HR  [  FT.  The  subroutines  also  set 
the  Initial  conditions  for  the  large-signal  equations. 

Subroutine  DIFFEQ  contains  Rowe's  differential  equations.  The 


computer  names  used  for  the  various  variables  are  given  In  Table  2. 
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the  equations  ire  store!  In  a  single  arrsv  called  FONF,  md  the  solu¬ 
tions  t>  the  equations  ire  stored  in  the  irrav,  VM.l'FS. 

Subroutine  DRIFT  perf  mi  the  same  tuties  for  DRFTF)  that  DIFF 
perform  for  DIFFFQ. 


C .  Integrals  In  Rowe’s  Differential  F.quatlons 

The  circuit  wave  amplitude,  the  circuit  wave  phase,  and  the  elec¬ 
tron  force  equations  each  contain  integrals  that  must  he  solved  at  every 
step  of  the  differential  equation  solution.  These  integrals  are  solved 
In  subroutines  FIST.  INTI.  AND  IST2. 


o  it  Mit.  t  Me  v  I N  r  solves  the  two  lntegrils  that  arise  from  the 
v  irt>-r  series  expansion  of  the  beam  charge  lensltv.  Since  the  electron 
Msxs  ire  evenlv  spiced  over  2*  radians  Initially,  the  Integration  over 
case  an  he  replinel  with  in  Integration  over  all  the  electron  disks. 
M  >re  'ver,  t  he  discrete  number  >f  electron  llsks  used  In  the  computer 
n  lei  ill  >ws  the  Integrals  to  he  replaced  bv  summations.  These  sumraa- 
■  :  >  is  ire  -inducted  nice  it  each  step  of  the  differential  equation 
s  1 :  ,»  i  m  . 

Sihriutfoe  1ST',  solves  the  Integral  from  the  space-charge  field 
expressl  m .  In  the  beam  charge  density  Integrals,  this  Integral  Is 

rep',  i  el  hv  a  summation  over  the  electron  disks.  The  space-charge 
sommatl  >n  must  he  conducted  for  each  electron  disk  at  each  step. 

sihrmtlne  INT2  solves  the  polynomial  approximation  to  the  space- 
•harge  Integral.  This  subroutine  runs  In  about  half  the  time  that  INTI 
takes  ti  run.  Since  the  space-charge  Integral  Is  evaluated  so  often, 
islng  the  polynomial  approximation  can  make  a  significant  difference  In 
run  time. 

0.  Numerical  Analysts 

Numerical  analysis  routines  are  used  to  solve  Rowe's  large-signal 
equations  and  to  calcualte  Bessel  functions  for  Pierce's  gain  parame¬ 
ter.  The  numerical  analysis  routines  are  DESOLV,  STEP,  INTRP,  MACHIN, 
and  BESSEL. 

The  differential  equation  solver  used  In  TWT1D  was  taken  from  a 
text  by  L.  P.  Shamplne  and  M.  K.  Gordon.®  The  solver  19  a  variable 
step,  variable  order  Adams  code.  It  consists  of  routines  DESOLV,  STEP, 


INTRP,  and  MACHIN.  Subroutine  DESOLV  Is  a  driver  that  calls  the  other 
subroutines  In  the  solver.  MACHIN  computes  the  machine  unit  roundoff 
error,  U,  which  Is  the  smallest  positive  number  such  that  1.0  +  U  > 
1.0.  Subroutine  STEP  Integrates  the  differential  equations  one  step, 
and  subroutine  INTRP  approximates  the  solution  at  the  end  point  by 
evaluating  the  polynomial  there.  The  code  Is  completely  explained  and 
documented  In  the  text  by  Shamplne  and  Gordon,  and  those  Interested  In 
the  code  should  consult  that  text. 

Subroutine  BESSEL  calculates  modified  Bessel  functions  of  the 
first  and  second  kind. 

E.  Program  Output 

As  Rowe's  equations  are  Integrated  down  the  traveling-wave  tube, 
various  parameters  describing  the  tube's  performance  are  printed  to 
output  data  and  plot  flies  and  displayed  on  the  terminal.  In  addition, 
variables  describing  the  beam  and  circuit  wave  must  be  printed  for  a 
restart  printout.  The  subroutines  that  handle  this  are  OUTPUT  and 
OUTRS. 

Subroutine  OUTPUT  controls  the  output  of  data  at  each  print  Inter¬ 
val.  It  prints  the  axial  position,  power,  and  phase  angle  to  the  plot 
file,  and  displays  the  position  and  power  on  the  terminal.  The  NPR 
variable  19  used  to  determine  which  variables  are  printed  to  the  output 
data  file.  In  addition,  OUTPUT  uses  the  power  and  lt9  derivative  with 
respect  to  position  to  calculate  the  saturated  power. 

Subroutine  OUTRS  controls  the  restart  printouts.  A  restart  output 
Is  used  when  an  engineer  wants  to  examine  the  effect  of  different  output 


-  29  - 


helices  without  having  to  run  the  Input  helix  In  each  case.  The  circuit 


wave  variables  printed  out  are  normalized  position,  Y,  circuit  wave 
amplitude,  A(n),  Its  derivative,  DA(n),  circuit  wave  phase,  THETA(n), 
and  Its  derivative,  DTHETA(n).  The  beam  variables  printed  out  are 
electron  phase  positions,  PHI(i,n),  and  electron  velocities,  U(i).  In 
addition,  If  the  polynomial  approximation  to  the  space-charge  function 
Is  used,  Its  coefficients,  CB(m)  are  printed  out.  These  data  contain 
the  Information  needed  to  restart  a  run  at  that  position. 

F.  Circuit  Changes 

Whenever  a  change  In  circuit  velocity  or  loss  occurs,  the  parame¬ 
ters  describing  the  circuit  must  be  recomputed.  Subroutine  ATTEN  han¬ 
dles  changes  In  loss  due  to  attenuators  or  severs,  and  subroutine  CHANGE 
handles  changes  in  velocity  due  to  velocity  steps  or  tapers. 

Attenuators  are  broken  into  20  sections  in  the  main  program.  In 
each  section,  subroutine  ATTEN  adjusts  Pierce's  loss  parameter,  calls  a 
subroutine  to  integrate  the  large-signal  equations,  and  checks  the  event 
code  list  to  see  If  any  velocity  changes  or  restart  printouts  occur  In 
the  section.  If  the  loss  In  a  section  exceeds  170  dB,  the  section  Is 
treated  as  a  sever,  and  the  drift  tube  equations  are  Integrated.  At  the 
end  of  the  sever,  the  subroutine  sets  Initial  conditions  and  Integrates 
the  large-signal  equations.  After  all  the  sections  are  integrated,  the 
loss  parameter  Is  reset  to  Its  original  value. 

Subroutine  CHANGE  first  checks  to  see  If  the  velocity  change 
Includes  a  change  In  helix  radius.  If  so,  the  fill  factor  is  rescaled, 
and  the  voltage,  beam  velocity,  and  plasma  frequency  are  recalculated. 


This  expression  is  plotted  along  with  Rowe  s  space-charge  function  for 


the  case  of  f(a'/b')  =  1.25  in  Fig.  9.  As  the  figure  shows,  the  agree¬ 
ment  between  the  two  functions  is  very  good.  Figure  10  compares  results 


Polynomial  approximation  to  Rowe's  space-charge  function 
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Fig.  10.  Power  plot  using  the  polynomial  approximation. 


of  a  TWT1D  run  using  the  full  space-charge  expression  to  a  run  on  the 
same  tube  using  the  polynomial  approximation.  The  results  of  the  two 
runs  are  almost  identical,  but  the  run  using  the  approximation  took  only 
half  the  time  of  the  run  using  Rowe's  full  space-charge  expression. 

Subroutine  SPACE  calculates  the  coefficients  for  the  polynomial 
approximation  using  the  function  SIMINT.  SIMINT  solves  the  integral  for 
the  coefficients  using  Simpson's  integration. 


•v-v-v-v  ■ 


V.  RESULTS  AND  CONCLUSION 


In  this  section,  the  validity  of  TWTID  is  investigated.  Two 
Teledyne  MEC  tubes  were  simulated  by  the  program,  and  the  results  are 
compared  to  the  measured  output  of  the  tubes. 

The  first  tube  modeled  is  an  I/J  band  pulse  tube.  It  has  6 
attenuators,  4  pitch  changes,  and  a  diameter  step.  The  tube  was 
simulated,  using  16  electron  disks  and  Rowe's  full  space-charge 
expression.  Several  runs  of  the  program  were  made  to  determine  the 
saturated  power  across  the  band.  The  results  of  the  program  and  the 
measured  output  power  are  plotted  in  Fig.  11.  As  the  figure  shows,  the 
agreement  across  the  band  is  very  good.  However,  at  the  low  end  of  the 
band,  the  predicted  power  is  1  to  2  dB  higher  than  the  measured 
output.  In  order  to  correct  this,  the  program  was  rerun  for  the  first 
three  data  points  with  the  harmonic  included.  Since  the  tube  is  a 
broadband  device,  harmonics  can  have  a  significant  effect  on  the  output 
power  at  the  low  end  of  the  band.  The  results  of  the  runs  with  the 
harmonic  included  are  shown  in  Fig.  12.  As  this  figure  shows,  the 
harmonic  does  lower  the  saturated  power  and  brings  the  2  curves  into 
closer  agreement  at  the  lower  frequencies. 
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The  second  tube  modeled  Is  an  E/H  band  CW  tube.  It  has  3 
attenuators  and  2  pitch  changes.  The  tube  was  simulated,  using  16 
electron  disks  and  the  polynomial  approximation  to  the  space-charge 
expression.  The  results  of  the  runs  and  the  measured  output  are  plotted 
in  Fig.  13.  Again,  the  agreement  between  the  two  curves  is  very  good. 

The  simulations  of  the  two  tubes  indicate  that  TWT1D  can  be  a 
useful  tool  for  the  design  engineer.  However,  care  must  be  taken 
whenever  a  nonlinear  program  is  used.  Many  of  the  parameters  that  are 
used  by  such  a  program,  such  as  the  fill  factor  and  the  axial  impedance, 
are  not  precisely  known  to  the  engineer.  Vet  these  parameters  have  a 
large  effect  on  the  output  of  the  program.  Figure  14  shows  the  results 
that  changing  the  fill  factor  can  have  on  the  output  of  the  program. 
Note  that  changing  the  fill  factor  from  0.4  to  0.6  moves  the  saturation 
position  as  much  as  2  inches. 

Perhaps  the  biggest  problem  with  any  large-signal  program  based  on 
Rowe's  approach  is  the  way  it  models  attenuators.  As  was  discussed  in 
Section  III,  attenuators  cause  oscillations  in  the  circuit  wave  phase 
angle.  This  problem  was  circumvented  in  TWT1D  by  replacing  attenuators 
with  severs,  but  that  is  not  a  complete  solution.  This  problem  should 
be  solved  before  more  elaborate  models  of  Rowe's  approach  are  attempted. 
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APPENDIX  A 


TWT  ID  SOURCE  CODE  LISTING 


FTN77  ,  S 
$CDS  ON 
i F I LES ( 3,3) 

$EMA  // 

SEMA  /TEMP/ 

PROGRAM  TWT 1 D < 3,91) 

C  LARGE  SIGNAL  ANALYSIS  PROGRAM  FOR  HELIX  TWT ’ S 

C  BASED  ON  THEORY  BY  J.E.  ROWE  S.  M.K.  SHERBA 

C  BY  THERESA  BRUNASSO 

C  TELEDYNE  MEC 

C  3165  PORTER  DRIVE 

C  PALO  ALTO,  CA  94304 

C  LAST  REVISED  < 860 5  1  5 . 07  1  7 > 

C  - 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

CHARACTER*64  I F I LE , 0 F  I  LE , P F I LE 

INTEGER  AT YPE! 6) , RSTART, SETU, SCS , SWITCH, SWLIST , TODAY!  15  )  , 
4TRIMLEN, TYPE 

DOUBLE  PRECISION  A  I S TEP ! 3 ) , A  I  TEMP ! 6 , 3 ) , CRTEMP f 6 )  , F I ! 3 )  ,  I COS  , 

&  I  FUNCT  ,  IS  IN,  10  ,  1 0GB  ,  I  1GB  ,  K1  ,  KL' ,  K3  ,  KOGB  ,  K1GB  ,  LA,  LAMAX  i  6  )  ,  LT  ,  MAG  , 
fcNORMZ , RFSTEP ( 3 ) , R  FTEMP  ! 6 , 3 ) , SHAPE ! 20 , 4 ) , TAPER ! 6 ) , VPS TEP ( 3 ) , 

& VPTEMP ( 6 , 3) , ZA( 19, 6) , ZATNE! 61 , ZATNS! 6) , ZC( 60 ) , Z I! 6) , ZLIST!  60  ■  , 
&ZRST( 10 ) 

COMMON  A(3),ANG(3),AXIMP(3),AXIMPN(60,3),A0!3),B!3),BETAf31, 
g>BETA0(3),BP,C(3),CB(5)  ,  CKTRAD  (  60  ),D(3),DA(3),DPC,DPDB,DTHETAi3,t 
fcD  F  f  3  )  , EBE AM ,EBLNC,EM,EPSI,ESAT.ETA, EWAVE  ,EXP1,F!31,FILL,G<3), 

&G AMA  f  3 )  , GS AT ,  ICOS ( 3 )  ,  I  FUNCT! 64 ) ,  INTSKP, ISIN!3),I0,K1,K2,K3, 

&  LA (20,6) , M AG ! 3 )  , N ATTN , NCHNG , ND , NE LECT , NORMZ , NPR , NUMF , 0 LD PD B , 0 LDZ , 
&PDB, PDB1 , PHI (64 , 3 )  , PHVEL! 3 ) , PH VE LN ( 60 , 3 ) , P I , P I N I ( 3 )  , PRNT , PS  AT , PW , 
&RFLSS ( 3 ) ,RFLSSN(60,3) , SCS , SETU , SKPINT, SMALL A , SMALLC , SWLIST ! 140) , 
&THETA ( 3 )  , TP  I , U ! 64 )  ,UTERM!64) ,UZER0,VP!3),VREL,V0,W(3)  .WP.WPWC, 
&YFINAL, YHALT, YPRINT, Z  ,  Z  F  ,  ZS  AT  ,  ZS  TOP  (1  40  ) 

COMMON  /TEMP/  FTEMP(64,64) 

DATA  SHAPE/  . 0 2 3 ,  . 0 3 3 ,  . 0 5 ,  . 067 ,  . 0 92 ,  .  1 2 ,  . 1 5 ,  .  1 93 ,  . 2 5 7 ,  . 3 3 3 ,  . 4 6  , 

&. 593,  .733,  .9,1.  ,.9,  .733,  .5, .2 33,  .033,  .013,  .030,  .066,  .129,  .231, 

& . 375 ,  .553,  .739,  .897,  .990,  .990,  .897,  .739,. 553,. 375,. 231,.  129,. 066, 
&. 03, .013, .033,  .233, .5, .733, .90, 1.  ,  .90,  .733,  .593,  .46,  .333,  .257, 
i. 193 ,  . 15 ,  . 12 , . 092, . 067, . 05,  . 033, . 02  3 , 1 . , 1 . , 1 .  , 1 .  ,  1 .  ,  1 .  ,  1 .  , 1 .  , 1  .  , 
fcl.  , I .  , 1 . ,  1.  ,  1  .  ,  I.  , I.  ,  1  .  ,  1.  ,  1 .  , 1 .  / 

CALL  FPARM! I F I LE , OF  I LE , P F I LE ) 

LU-  1 

PI =4 . FATAN (  1  .  ) 

EM  -  1 . 7588047D  1  1 
EPS  I  =  8 . 85418782D-  1 2 
Kl-.5/(PI*EPSI*SQRT(2.*EM) ) 

K2= . 5/EMF2 . 997924 5800 8* »2 
K3- . 5/K2 

SMA L LC -  2 . 9979245RD08F39.37 

NCHNG  -0 
NATTN -  0 
SETU  0 
FIRST  0 
ND  =0 

DO  1-1,64 

FTEMP i  I  ,  I  )  =  0 
END  DO 


’  EM-e/m  !Coul/kg)' 

'  PERM  ITIVI  TTY  OF  FREE  SPACE  t'F/ni 
!Kl=l/(2*PI*EPSI*SQRT!2*e/m) ) ' 

' K2  =  ( m*c**2 ) /  (  2*e )  ' 

'SPEED  OF  LIGHT  lin/s)1 


-  41  - 


'  s'  *.*  ‘ 

.A  ’  *.  *  ■ 


A  •.*, 


OPENING  DISPLAY 


CALL  F  T  I M  E  ( TODAY ) 

WRI  TE<  LU  ,  10  )  (  TOD  AY  (  I)  ,  I  i  ,  1  5  ) 

10  FORMAT;  ,/  15X TWT1D  ’  ,  1  5  A  2  ,  ' 

IF  (IF!  I.K  (  1  :  1  >  .  FO  .  CHAR  O')  THEN 
WR  ITE (  LU , 20 ) 

20  FORMATf  //’  INPUT  DATA  FILE  NAME’'  ') 

READ! LU, 30)1  FILE 
30  FORMAT; A) 

END  IF 

IF  (OFILEM:  1)  .EQ.  CUARiO't  THEN 
WRITE ( LU , 40 ) 

40  FORMAT;//’  OUTPUT  FILE  NAME0  _’) 

READ  f LU , 30 ) OF  I LE 
END  IF 

IF  i  P  F  I  L  E  <  I  :  1  j  .EQ.  CHARiO7))  THEN 

write; LU, 50 ) 

50  FORMAT i / / ’  POWER  PLOT  FILE  NAME0  ’) 

RE  AD i  LU, 30  1PFILE 
END  If 

WRITE;  LU  ,  60  :  OF  I  LE  (  I  :  T  R  I  M  I.  E  N  (  0  F  I  LE  )  i  ,  PF  I  LE  1  :  TR  IMI.EN  PF  I  LE 
60  FORMAT;/  ’  ALL  OUTPUT  DATA  WILL  BE  ROUTED  TO  FILE  ’,A 
f.'  ALL  PLOT  DATA  WILL  BE  ROUTED  TO  FILE  ’,A) 


READ  DATA  FROM  INPUT  FI  IE 


0  P  E  N  <  9  9  ,  F  I  L  E  I  F  I  L  E  > 

READ! 99, *  , V 0,  10, RAD, FILL, TP  I, ZF 
READ  <  99, *  NPMF, NUMC , NI  MA , NUMR 
DO  1=1, NUMF 

READ ( 99, *  i  FI '  I  -  , P INI (  I  )  , P  H  V  E  L  <  I  .  , AX  IMP (  I  ) , RFLSS (  I)  , THETA.  I 
END  DO 

IF  (NUMC  .GT.  0)  THEN 
DO  1=1, NUMC 

READ(99,*)7, 1(  I  )  ,  ORTEMP  (  I), TAPER)  I) 

END  DO 

RE AD ( 99 , *  )  ( (VPTEMP; I , J ) , 1=1, NUMC ) , J= 1 , NUMF) 

READ; 99, * ) ( ( A  I  TEMP;  I,J;  ,  1  =  1, NUMC ) , J= 1 , NUMF ) 

RE AD ( 99 , *  )  (  iRFTEMP;  I , J )  ,  1  =  1, NUMC )  , J  =  1 , NUMF ) 

END  IF 

IF  ( NUMA  .GT.  0 !  THEN 
DO  1=1, NUMA 

READf  99, * 'ZATNSf  I  i , ZATNEII) , LAMAX( J ) , ATYPEf  I) 

END  DO 

DO  1=1, NUMA  ’SET  UP  ATTENUATOR  SECTIONS' 

ASTEP=(ZATNE;I)-ZATNS(I))*.05 
ZA(  1  ,  I) =Z ATNS I  I  )  +  ASTEP 

la; l , i ) =lamax; i ) *shape; l , atype; n ) 

DO  K  =  2  ,  19 

la;k,i)=lamax;i)*shape(k,atype;i)) 

ZA;k, I )  -  Z  A i  K-  1  ,  I ) *ASTEP 
END  DO 

la; 20,  i )  =  lamax;  i )  asiiape;  20 ,  atype;  i ) ) 

END  DO 
END  IF 
NUMT  =  0 
K  =  0 

IF  (NUMC  .GT.  0)  THEN 


•SET  UP  TAPERS' 


DO  I  -  1  ,  N'UMC 
1  M  1  I  -  1 

IF  I A  P  F H  (  I  1  . EQ .  0)  THEN  'PITCH  CHANGE* 
k  -  K  -  1 

NC  N  U  M  T  *  1  0  +  K 
X  C  (  N  C  i  -  Z  I(  I  ) 

CKTHAD ( NC 1  C  R  T  E  M  P (  I ) 

DO  J  -  1  , NUMF 

PHVELNt  NC , J I  -  V  PTEMP  (  1  ,  .1) 

A  X  I MPN  (  NC , J ) =  A  I  TEMP  (  I  ,  .7  > 

RFLSSNI NC , J) =RFTEMP( I , J ) 

END  DO 

ELSE  'VELOCITY  TAPER’ 

NUMT  =  NUMT-t-  1 
TSTEP-TAPERf  H  *  .  1 
N’T  =  I  +  1  0*  (  NUMT- 1  )  *K 

zc : NT ) - Z I : I) 

IF  (I  .EQ.  11  THEN  'FIRST  PITCH  CHANGE' 

DO  J= 1 , NUMF 

VPSTEPI J )  =  ( V PTEMP (  1  ,  J  1  -PHVELI J )  ) *  .  1 
A  I  STEP! J )  =  ( A  I TEMP!  1  ,  J  )  -  AX  IMP ( J )  I  *  .  1 
RFSTEPI J ) - ( RFTEMP ( 1 . J ) -RFLSSI J 1) * . 1 
PH VELN( NT, J ) -PHVELI J ) + VPSTEPI  J i 
AXIMPN(NT,J ) =  AX  IMP ( J )  +  A  I  STEP (  J  1 
RFLSSN(NT,J1 =  RF  LSS ( J ) +  R  FS  TE  P (  J  ' 

END  DO 
ELSE 

DO  J=  1  , NUMF 

VPSTEPI  J  )  =  (  VPTEMP  (  I  ,  J  1  -  VPTEMP  <  IM1  ,  .1  :  *  .  1 
A  I  STEP! J )  =  ( A  I TEMP ( I , J ) -AITEMPf IM1  ,  J )  ) * .  1 
RFSTEPI J)=( RFTEMP ( I  ,  J  ) -RFTEMP  < IM1  ,J)  )*.  1 
PH VELN( NT, J )  = VPTEMP! IM1 , J 1 + VPSTEP ( J ) 

AX  I MPN ( NT, J )  =AITEMP( IM 1 , J )  +  A  I STEP ( J) 
RFLSSNI NT, J ) = RFTEMP ( IM1 , J) +RFSTEP  J) 

END  DO 
END  IF 

CKTRAD ( NT ) =  C  RTEMP ( I ) 

DO  NSTEP=2 , 10 
DO  J=  I  , NUMF 

NT= 1 0  * ( NUMT -  1 ) +NSTEP+K 
NTM 1  -  NT-  1 

ZCl NT ) -Z I( I )+( NSTEP-1 ) *TSTEP 
CKTRAD ( NT) =CRTEMPI  I  ) 

PHVELN(NT,J)=PHVELN(NTM1  ,J  1+VPSTEP(.J  ) 

AX IMPN ( NT , J ) -AX IMPN ( NTM1 , J 1 + A  I  STEP ( J  t 
RFLSSNf  NT,  J)=RFLSSN(NTM1  ,  J  )  ♦RFSTKPi'  .1  1 
END  DO 
END  DO 
END  I F 
END  DO 
END  IF 

IF  ( NUMR  .GT.  0)  RE  AD ( 99 , * ) ( Z RST ( II  ,  I ^ 1 , NUMR  1 
READ ( 99, * ) SCS , NELECT, DPC, NPR, INTSKP, RSTART, F DUMMY 
IF  (RSTART  .EQ.  0)  CLOSE(99) 


COMPUTE  ROWE  VARIABLES  AND  PIERCE  PARAMETERS 


V2 -  VO* ( 1 . - I0*V0**(-1 . 5 ) *K 1  * ( . 5-LOGC FILL )  )  ) 

VREL-K2*(  1  .  1  .  /  (  1  .  «-K.1*V2)**2)  'RELATIVISTIC  VOLTAGE' 


i  m  i  nj '  w* 


■  'iwn  mm  ww^twi 


fc 


>. 

v; 


R  I  =  F  LO AT (  I  ) 

PHI ( I , N ) =2 . *RN*P I *<  RI- I  .  )  /FLOAT  (  NELECT)  +  (  KN  *  THETA  (  1  ) 

&  -THF.TA'N)} 

END  DO 
END  DO 
ELSE 

DO  N=  1  ,  NUMF 

M=  F  I  ( N ) / FDUMMY 
DO  1=1, NELECT 

PHI ( I , N) =2*M*PI * ! I-l>  'NELECT -THETA  (  N) 

END  DO 
END  DO 
END  IF 

DO  N= 1 , NUMF  'INITIAL  COND .  FOR  CIRCUIT  WAVE! 

G(NI=-D(N)*A(N)*( l.+C<N)*B(N) )*W(N)*C(N)/WlCl 
F<  N ) =-B ( N) *W( N) *C(N1 /W1C 1 
END  DO 
END  IF 

PW=2.*C!l)*I0*VREL*A(l)*A!l)*(l.-C(l)*F(l))/<l+C(l)*B(l)) 

OLDPDB= 1  0  *  LOG 10 ( PW ) +30 

0LDZ=Y/N0RMZ 

SKP I  NT  =  0 

C  PRINT  OUT  VALUES  CALCULATED  FROM  INPUT 

C  - 

OPEN! 98, FILE=OFILE  ) 

OPEN ( 97 , F I LE  =  PF I LE ) 

IF  ( NPR  -EQ.  4)  0PEN(96, FILE= *  ELECTRON. PLT*  ) 

WRITE! 97, 3 )OLDZ, (PINI(N),N=1, NUMF ) , < THETA ( M) , M=1 , NUMF ) 

3  FORMAT! 1X,F7.4,6(2X, F8 . 4  ) ) 

CALL  FTIME ! TODAY ) 

WRITE (98, 101 ( TODAY ( I ) ,1=1,15) 

WRITE (98, 5) IFILE( 1 : TRIMLEN ( I K I IE ) ) 

5  FORMAT(20X, ’ RESULTS  FOR  INPUT  FILE  ’,A//) 

WR I TE ( 98 , 15) 

15  FORMAT ( ’  N  FREQ.(gHz)  INPUT  PWR.(dBm)  REL.  PHASE(rad)’) 

DO  N  = 1 , NUMF 

WRITE (98, 25) N, FI (N) , PIN  I (N) , THETA (N) 

25  FORMAT (  IX,  I1,3X,F8.4,7X,F8.4,  10X.F8.4) 

END  DO 

WRITE (98, 35)  10 , VO , VREL, UZERO/SMALLC , TP  I , WP , R 
35  FORMAT(/15X,  ’CATHODE  CURRENT ( Amps ) =  ’  , F7 . 4 / 1 6X ,  ’ HE L I  X  VOLTAGE! Volt 
&s)=  ’, F9. 2/ 12X,  ’ EFFECTIVE  VOLTAGE ( Vo  1 1 s ) =  ’  , F9 . 2 / 1 2 X ,  ’ NORM A L I Z ED  D 
&EAM  VELOCITY=  ’, F7 . 4/20X ,’ HELIX  PITCH(tpi)=  ’, F6 . 3/ 12X ,’ PLASMA  FRE 
&QUENCY ! rads/s ) =  ’, E 1 1 . 4/3X PLASMA  FREQUENCY  REDUCTION  FACTOR^  ’ ,F 
8.7  -  4 ) 

WRITE(98,45)QC,SMALLA,FILI.  ,-DPC,  NE  LECT  ,  NUMA  ,  NUMC 
45  FORMAT( 14X,  ’ SPACE  CHARGE  PARAMETER=  ’  , F 7 . 4/ I  5 X ,  ’ ME AN  HELIX  RADIUS! 
&in)=  ’  , F8.5/20X,  ’ BEAM  FILL  FACT0R=  ’  , F7 . 4 / 1 0 X ,  ’CIRCUIT  PRINT  INTER 
AVAL(  in ) =  ', F7.4/12X,  ’NUMBER  OF  ELECTRON  DISKS=  ’,  12/ 15X ,’ NUMBER  OF 
&  ATTENUATORS=  ’, 12/ 19X ,’ NUMBER  OF  CHANGES=  ’,12/) 

WRITE ( 98 , 55 ) 

55  FORMAT  (  ’  N  NORM.  PHASE  VEL.(Vp./c)  AXIAL  IMP.  (ohms)  GAIN  PARAMET 

&ER  ’  ) 

DO  N= 1 , NUMF 

WRITE! 98, 65 )N, VP ( N ) /SMA LLC , AX  IMP ( N ) , C ( N ) 

65  FORMAT ( IX, I1,6X,F8.4, 13X.F8.4, 13X.F5.4) 

END  DO 


V.'- 


y.  • 

v-/. 

fl. 

Orr; 


-  45  - 


o  o  o  o  u 


SORT  ALL  EVENTS  INTO  AXIAL  POSITION  LIST 


DO  1=1,140 
ZSTOP( I ) =0 
SWLIST( I)=l 
END  DO 
NSTOP=0 

IF  ( NUMR  .GT.  0)  THEN  !SORT  RESTART  DATA  POSITIONS! 

DO  1=1, NUMR 

ZLIST( I ) =ZRST( I ) 

END  DO 

NLIST=NUMR 

TYPE=2 

CALL  ZSORT( TYPE, Z LIST, NLIST, ZSTOP, NSTOP, SWLIST) 

END  IF 

IF  (NUMA  .GT.  0)  THEN  !SORT  ATTENUATOR  START  POSITIONS 

DO  1=1, NUMA 

ZLIST( I ) =ZATNS( I ) 

END  DO 

NLIST=NUMA 

TYPE=3 

CALL  ZSORT ( TYPE, Z LI  ST, NLIST, ZSTOP, NSTOP, SWLIST) 

DO  1=1, NUMA  ! SORT  ATTENUATOR  STOP  POSITIONS! 

ZLIST( I ) =Z ATNE ( I ) 

END  DO 
TYPE=5 

CALL  ZSORT ( TYPE, Z LIST, NLIST, ZSTOP, NSTOP, SWLIST) 

DO  1=1, NUMA  ! SORT  ATTENUATOR  SECTIONS! 

DO  J=l,  19 

ZLIST( J) =ZA(J, I) 

END  DO 
NLIST= 19 
TYPE  =  1 

CALL  ZSORT ( TYPE, ZL I  ST, NLIST, ZSTOP, NSTOP, SWLIST) 

END  DO 
END  IF 

IF  ( NUMC  .GT.  0)  THEN  !SORT  CIRCUIT  CHANGES! 

NLIST=NUMC+9*NUMT 
DO  1=1, NLIST 

ZLIST(  I  ) =ZC( I) 

END  DO 
TYPE=7 

CALL  ZSORT ( TYPE, Z LIST, NLIST, ZSTOP, NSTOP, SWLIST) 

END  IF 

IF  ( RSTART  .EQ.  1)  THEN  ! FI ND  RESTART  POSITION' 

DO  N=l, NSTOP 

SWITCH=SWLIST(N) 

CALL  EVENT (SWITCH, ICHNG , IRSTRT , IASTRT, IASTOP) 

IF  (IRSTRT  .EQ.  1)  ND  =  N 
END  DO 
END  IF 


PROGRAM  CONTROL  CODE 


DO  WHILE  ( ND  . LE .  NSTOP) 
ND=ND+1 

YHALT=ZSTOP(ND)*NORMZ 
CALL  INTSUB(Y, FIRST) 
SWITCH=SWLIST( ND ) 


CALL  EVENT (SWITCH, ICHNG , IRSTRT , IASTRT, IASTOP) 
IF  (IRSTRT  .EQ.  1)  CALL  OUTRS(Y) 

IF  (ICHNG  .EQ.  1)  CALL  CHANGE ( Y ) 

IF  (IASTRT  .EQ.  1)  CALL  ATTEN( Y , FIRST, IASTOP ) 
END  DO 

YHALT=ZF*NORMZ 
CALL  INTSUB(Y, FIRST) 

SSGAIN^PDBl-PINI (  1  ) 


ENDING  SUMMARY 

IF  ( PDB 1  .GT.  PSAT)  THEN  ! TUBE  HAS  NOT  SATURATED! 

PS AT=PDB I 
GSAT=PDB1-PINI ( 1) 

ES AT-ETA* 100. 

ZS AT=Z 
END  IF 

WRITE(98,70)PSAT,GSAT,ESAT,ZSAT,  SSGAIN 
70  FORMAT(/26X, ’ENDING  PRINT  OUT ’/ 1 9 X ,* SATURATED  POWER  (dBm)=  ’.F7.3/ 
&17X, ’GAIN  AT  SATURATION  (dB)  =  ’, F7 . 3/ 15X ,’ SATURATION  EFFICIENCY  ( X 
&)=  ’, F7 . 3/I5X, ’ LENGTH  AT  SATURATION  (in)=  ’ , F7 . 3/ 18X , ’ SMALL  SIGNAL 
&  GAIN  ( dB ) =  ’ , F7 . 3 ) 

CALL  FT IME ( TODAY ) 

WRITE (98, 10 ) ( TODAY ( I ) , 1= 1 , 15) 

C  LOSE ( 98 ) 

C  LOSE ( 97 ) 


SUBROUTINE  SPACE(AB,BB,SCS,R,EXP1,CB,PI) 

COMPUTES  THE  SPACE  CHARGE  REDUCTION  FACTOR  (R),  THE 
DERIVATIVE  OF  THE  LOG(l-R),  AND  THE  COEFFICIENTS  (CB)  FOR 
THE  POLYNOMIAL  APPROXIMATION  TO  THE  SPACE  CHARGE  FUNCTION 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS 

DOUBLE  PRECISION  I0AB,I1AB,K0AB,K1AB,I0BB,I1BB,K0BB,K1BB,CA(5), 
&CB ( 5 ) 

EMA  SCS.EXPl.CB.PI 

CALL  BESSEL(I0AB, 1 1 AB , KOAB , KlAB , AB ) 

CALL  BESSEL( IOBB,  I1BB,K0BB,K1BB,BB ) 
R0=SQRT(ABS(1-BB*(K0AB*I1BB+I0AB*K1BB)/I0AB)) 

FO  =  LOG ( 1-RO) 

DELTA= . 001 
BB=BB+DELTA 

CALL  BESSEL(  IOBB,  I1BB,K0BB,K1BB,BB) 

R=SQRT ( ABS ( 1-BB*(K0AB*I1BB+I0AB*K1BB)/I0AB)) 

FI -  LOG ( 1-R) 

BB=BB+DELTA 

CALL  BESSEL( IOBB, I1BB,K0BB,K1BB,  BB) 

R=SQRT( ABS( 1 -BB * ( KO AB* I 1BB+ 10 AB*K 1 BB ) / 10 AB ) ) 

F2  =  LOG ( 1-R ) 

F--(  -  F2  +  4*F1-3*F0) / (2*DKLTA)  'CENTRAL  DIFFERENCE  FORMULA 

R  =  RO 

BB-BB-2*DELTA 

EXP1-EXP(BB*F) 

IF  (SCS  .EQ.  2)  THEN  ! PO  LYNOM I A  L  APPROXIMATION' 

DO  1-1.5 

C  A  (  I  )  =  0 


APPROXIMATION' 


o  n 


CB { I ) =0 
END  DO 
DO  1=1,3 
N  =  2 *  I- 1 

C  A  ( I)=2. /PI*SIMINTIN, EXP1) 

C  A ( I) =CA( I) -SORT ( .002) /PI 
END  DO 

CB (  1  ) =CA( 1  )-3. *CA(2) +5 . *C A ( 3 ) - 7 . *C A ( 4 ) +9 . *C A ( 5 ) 
CB(2)=4. *CA(2)-20. *CA(3)+56. *CA(4)-120. *CA(5) 

CB ( 3) =16. *CA(3)-112. *CA(4)+432. *C A ( 5  ) 

CB(4)=64. *CA(4)-576. *C A ( 5  ) 

CB(5)=256. *CA(5) 

END  IF 
RETURN 


FUNCTION  SIMINT(N, EXP1 ) 

USES  SIMPSON  INTEGRATION  TO  COMPUTE  THE  COEFFICIENTS  FOR  THE 
POLYNOMIAL  APPROXIMATION  TO  THE  SPACE  CHARGE  FUNCTION 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  LOWLIM 
EMA  EXP1 

SPCHI Y.EXPl )=.5/PI*( . 5*(PI-Y)-2. * ATAN ( S IN ( Y ) / ( EXP1-COS ( Y ) ) ) 

8,  +ATAN(SIN( Y)/(EXP1**2-C0S(Y) ) ) )  !  ROWE  ’  S  SPACE  CHARGE  FUNCTION! 

CHEV(N,X)=COS(N*ACOS(X) )  ! CHEB YSHEV  POLYNOMIALS! 

F(X)=SPCH(PI*(X  +  1.  ) , E  XP 1 ) *CHEV(N,X1/SQRT( 1 . -X*X)  ! INTEGRAND ! 

P I =4 . *ATAN (  1  .  ) 

DE  LTA= 1 . D-04 
LOWLIM=-l . +DELTA 
UPL IM= 1 . -DELTA 
S  UM2  =  0 
SUM4  =  0 

X=LOWLIM+DELTA 
DO  WHILE  (X  . LT.  UPLIM) 

SUM4=SUM4+F ( X ) 

SUM2=SUM2+F ( X+DELTA) 

X  =  X  +  2 . *DELTA 
END  DO 

SIMINT=DELTA/3. *(4. *SUM4+2. *SUM2+F ( LOWLIM ) +F ( UPLIM ) 

8,  +4.  *F(UPLIM-DELTA)  ) 

RETURN 

END 


SUBROUTINE  BESSEL ( A , B , C , D , X ) 

RETURNS  MODIFIED  BESSEL  FUNCTIONS  OF  1ST  AND  2ND  KIND 
A=  10 ,  B  =  1 1 .  C  =  K0 ,  D  =  K 1 ,  X  =  ARGUMENT 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 
IF  IX  . LT.  1 . 0D-10  )  X=1 . 0D-10 


. 5772157  +  LOG ( X* . 50 ) 
:X*X*. 25 


o  o  o  o  o  o 


10  G4=G4+1 . 

G3=G3+ 1 . /G4 
G1=G1*G6/(G4*G4 ) 

G2=G2*G6/(G4*(G4+1 . ) ) 

A  =  A+G 1 
B  =  B  +  G2 
C  =C+G 1 *G3 

IF  ( G 1 *G3  .GT.  1.0D-06)  GO  TO  10 

C=C-A*G5 

B= . 5*B*X 

IF  (X  .GE.  2. )  THEN 
X 1 =2 . /X 

CK=X1*5. 3208D-04  -  2.5154D-03 
CK=Xl*CK+5.87872D-03 
CK=X1*CK-1 . 062446D-02 
CK=Xl*CK+2. 189568D-02 
CK=Xl*CK-7. 832358D-02 
CK  =  X 1 *CK+ 1 . 25331414 
C=CK/(SQRT(X)*EXP(X) ) 

END  IF 

D  = ( 1 . /X-B*C)/A 

RETURN 

END 

C=-  =  =  --  =  =  =  =  =  =  =  =  =  -  =  --  =  =  =  =  -  =  =  =  =  -  =  =  -=  =  =  ^-  =  =  =-z:=.--  =  =  =  z:  =  zz--  =  =  -  =  --  -  = 

$EMA  // 

SUBROUTINE  I  NTS UB ( Y , F I RST ) 

INSURES  THAT  THE  APPROPIATE  PRINTOUTS  ARE  CONDUCTED, 

AND  THAT  THE  INTEGRATION  GOES'  TO  THE  END  OF  THE  INTERVAL 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS , SWITCH, SWL 1ST , SETU 

DOUBLE  PRECISION  ICOS , I FUNCT , IS  IN , 1 0 , K 1 , K2 , K3 , LA , MAG , NORMZ 
COMMON  A ( 3 ) , ANG ( 3 ) , AXIMP(3) , AXIMPN(60, 3) , A0(3) , B(3) , BETA (3) , 
&BETA0 ( 3 ) ,BP,C(3),CB(5) ,CKTRAD(60) ,D(3),DA(3), DPC  ,  DPDB , DTHETA ( 3 ) , 
&DF(3) , EBEAM.EBLNC, EM, EPS  I , ES AT , ETA , EWAVE , EXP1 ,F(3),FILL,Gf3), 
&GAMA ( 3 ) , GS AT , ICOS(3) , I FUNCT (64), I NTSKP , I S I N ( 3 ) , I 0 , K 1 , K2 , K3 , 

&LA( 20,6) , MAG ( 3 ) , N ATTN , NCHNG , ND , NELECT , NORMZ , NPR , NUMF , OLD PD B , 0 LDZ , 
&PDB,PDB1,PHI(64,3) ,PHVEL(3) , PH VE LN ( 60 , 3 ) ,PI,PINI(3) , PRNT , PS AT , PW , 
&RFLSS(3) ,RFLSSN(60,3) , SCS , S ETU , SKP I  NT , SMALL A , SMALLC , SWL I  ST (  140 )  , 
&THETA ( 3 ) , TP I, U 164) ,UTERM(64) , UZERO , VP ( 3 ) , VRE L , VO , W ( 3 ) .WP.WPWC, 
AYFINAL, YHALT, YPRINT , Z , ZF , ZSAT, ZSTOP( 140) 

IF  (YHALT  .GT.  YFINAL)  YHALT=YFINAL 
DO  YPRINT=Y+PRNT, YHALT, PRNT 
CALL  DIFF(Y) 

Z=Y/NORMZ 

CALL  OUTPUT( FIRST) 

END  DO 

IF  (Y  .LT.  YHALT)  THEN 
YPRINT=YHALT 
CALL  D I FF  C  Y) 

Z= Y/NORMZ 

CALL  OUTPUT( FIRST) 

END  IF 
RETURN 
END 


SUBROUTINE  ZSORT(TYPE,ZLIST,NLIST,ZSTOP,NSTOP,SWLIST1 
SORTS  THE  POSITIONS  AT  WHICH  CIRCUIT  CHANGES,  ATTENUATORS, 
ETC.  OCCUR,  AND  FORMS  THEM  INTO  A  SINGLE  LIST. 
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ZLIST=LIST  OF  STOP  POSITIONS,  NLIST=LENGTH  OF  ZLIST, 

ZSTOP  =  SORTED  LIST  OF  STOP  POSITIONS,  NSTOP=  LENGTH  OF  ZSTOP 
SWL I ST=  L 1ST  OF  EVENT  TYPES  CORRESPONDING  TO  ZSTOP 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SWLIST, TYPE 

DIMENSION  ZLIST(NLIST) ,ZSTOP( 140) ,SWLIST( 140) 
EMA  ZSTOP, SWLIST 
DO  I = 1 , NLIST 
N=  1 


DO 

WHILE 

!  ZLI 

ST  (  I ) 

.GT. 

ZSTOP ( N ) ) 

N  =  N+- 1 

IF  (N 

.GT. 

NSTOP)  GO 

TO 

200 

END 

DO 

IF 

(ZLIS 

T(  I) 

.  NE. 

ZSTOP 

(N)  ) 

THEN 

DO  J  = 

NSTOP 

,  N ,  -  1 

ZSTOP! J 

+  1 ) =ZSTOP  f 

J) 

SW 

LIST! 

J  +  I  )  = 

SWLIST! J) 

END  DO 
SWLIST(N) =1 
END  IF 

200  IF  ( NSTOP  .EQ.  0)  N=1 
ZSTOP(N)=ZLIST(  I) 

NSTOP  =  NSTOP-*- 1 
SWLIST(N)=SWLIST(N)*TYPE 
END  DO 
RETURN 
END 
C 

SUBROUTINE  E V E NT ! S W I TC H ,  ICHNG ,  IRSTRT,  IASTRT,  IASTOP) 

C  CHECKS  SWLIST  TO  DETERMINE  WHICH  EVENTS  OCCUR. 

C  EVENT  TYPE  CODE:  2--PRINT  RESTART  DATA;  3--ATTENUATOR  START; 

C  5--ATTENUATOR  STOP;  7--CIRCUIT  CHANGE 

C  - 

INTEGER  SWITCH 
ICHNG -  0 
I RSTRT  =  0 
I ASTRT^O 
IASTOP=0 


TEST=MOD ! 

SWITCH, 

2) 

IF  (TEST 

•EQ.  0) 

I RSTRT= 

TEST=MOD ( 

SWITCH, 

3) 

IF  ! TEST 

.EQ.  0) 

I ASTRT= 

TEST=MOD ( 

SWITCH, 

5) 

IF  ! TEST 

.EQ.  0) 

I ASTOP- 

TEST=MOD ( 

SWITCH , 

7) 

IF  (TEST 

.EQ.  0) 

ICHNG= 1 

RETURN 

END 


$  EM A  // 

SUBROUTINE  CHANGE(Y) 

C  RECOMPUTES  CIRCUIT  PARAMETERS  WHEN  A  CHANGE  IS  REACHED 

C  - 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS, SWITCH, SWLIST, SETU 

DOUBLE  PRECISION  ICOS.IFIJNCT,  ISIN,I0,KI,K2,K3,I0GB,I1GB,K0GB,K1GB, 
ALA, LT , MAG , NORMZ 

COMMON  A  (  3  I  , ANG ( 3 ) , AX  IMP! 3 ) , AX IMPN ( 60 , 3 ) , AO ( 3 ) , B ( 3 ) , BETA ( 3  )  , 


fUSETAO  (3)  ,  BP,C(3)  ,CB(5  )  ,  CKTHAD  (  60  )  ,  IM  3  '  ,  D  A  (  3  )  ,  DPC  ,  DPDR  ,  D  THETA  <  3  )  , 
&DF(  A  )  ,  EBEAM,  EBLNC,  EM,  t'PSI  ,  E  S  A  T  ,  E  T  A  ,  EWAVE.EXPl  ,  F  (  3  )  ,  FI  LL,  G(  3  )  , 
K.GAMA(  3  )  ,  GSAT  ,  I  COS  (  A  )  ,  I  FCNCT  t  64  )  ,  INTSKP  ,  IS  INI  3  )  ,  10  ,  K  1  ,  K2  ,  K3  , 

*<LA(  20 , 6  )  ,  MAG(  3  )  ,  NATTN  ,  NCHNG  ,  ND  ,  SELECT,  NORMZ  ,  NPR  ,  NUMF  ,  0  LD  PDB  ,  C  Ml 
R.PDB  ,  PDB1  ,  PHI  i  64 , 3  )  ,  PHVEL(  A  '  ,  PHVELN.  60 , 3  )  ,  P  I  ,  PINI  (  3  1  ,  PUNT  ,  PS  AT  ,  P 
8.RFLSS  i  3  )  ,RFLSSN(60,3)  ,  SCS  ,  S  ETl: ,  SKP  I  NT  ,  SMALL  A  ,  SMALLC  ,  SW  L  I  ST  (  140  ,  , 
A  THETA  (  3  )  ,  TP  I  ,  U  <  64  1  ,  UTERM  (  64  I  ,  11  ZERO  ,  VP  ,  A  1  ,  VREL  ,  VO  ,  W(  3  )  ,  WP  ,  WPWC  , 
&YFINAL,YHALT,YPRINT,Z,ZF,ZSAT,ZSTOP( 140' 

NCHNG-NCHNG+ 1 
ANEW^CKTRAD ( NCHNG  ) 

IF  (ANEW  . NE .  SMALLA)  THEN  'NEW  HELIX  RADIUS! 

F r LL- F I LL*SMALLA/AN£W  'FILL  FACTOR' 

V2  =  V0*(  1  . -I0*V0**( -1 . 5  )  *K1*(  .  5-I.OGi  FI  LL)  )  ) 


PHASE  VELOCITY' 


VRELrK2*( !.-!./(  l.+K3*V2)**2)  'RELATIVISTIC  VOLTAGE 

UZERO  -  SORT ( VREL ) *2 . 33502D  +  07  'BEAM  VELOCITY  dn/si 

DENOM-PI*EPSI*  UZERO*  B  P*B  P 

WP -  SORT i EM*I0*61024. /DENOM)  'PLASMA  FREQ,  (rads/s' 

END  IF 
DO  J  -  1  , NUMF 

IF  (PHVELN (NCHNG, 1 )  . LT.  1)  THEN  'PHASE  VELOC 

VPNEW=PHVELN( NCHNG, J ) *  S  M  A  L  L  C 
RETAN  U-VP(  J)  *  B  E  T  A  (  .1  )  /VPNEW 
GAMANU-SQRT(  BETANU**2-8ETA0i  .1  i  **2  > 

ELSE  '  T  P I  ' 

GAMANU=GAMA ( J ) * ANEW *P H VE LN ( NCHNG  ,  1  )  / ( SMALL A*TP I 
BET  ANII- SORT  (  GAMANU«*2  *  BETAO  (  I  >*«2( 

V  P  N  E  W  =  W ( J ) /BETANU 
TPI  - P H V E L N ( NCHNG,  1  i 


END  IF 

Zl-AXIMP(J)  'OLD  CKT  IMPEDANCE' 

IF  ( AX  IMPNf NCHNG , J  )  .EQ.  99)  THEN  'TPI  SCALING' 

AX  IMP(  J  )  =  AXIMP<  .1  )  *(  BETA(  J  i  /RKTANU)  *«3*(  GAMANU/GAMA(  J  )  '  **4* 
(  EXP  I  6664*GAMANU*ANEW  !  /  EXP  .  6664*GAMA(  .1  )*SMALLA)  )**3 
ELSE 


AX’  IMP(  J  )  --  AX 
END  IF 
Z  2  ■  A  X  I M  P  <  J  ) 
SMALLA-ANEW 
GAMA  t J I  - G AMANU 
RETA  (  J  )  --  BETANU 
VP  r  J  )  - VPNEW 


AX  IMPS' NCHNG , J ) 


•NEW  CKT  IMPEDANCE' 

' HELIX  RADIUS  (in'' 

' RADIAL  PROP.  CONST.  • 

’  PROPAGATION  CONST  .  1 
' PHASE  VELOCITY  (  in  s' 


CALL  BESSEL  I  10 GB,  I  1GB, KOGB, K 1 G  B , G  A  M  A ( J) *BP) 


ZAVG-AXIMP' J  )*(  I0GB**2-  I1GB**2 
C  1  J  )  -  I  I  0  *  Z  A  V  G  *  .  2  5  /  V  R  E  L  *  *  (  1  .  .3 
B  (  .1  >  -  (  UZERO/  VP(  J  )  -  1  .  i  ,  C  ( .1  ) 
RFLSS  I  J  )  -RFLSSN  NCHNG  .  .1  i 
LT  RFLSS ( J  )  / LOG  1  0 1  KXP<  1  .  )  ) 

D ( J ) - LT* . 05/ ( W( J  )  /  UZERO *C ■  J 
A(J)-A(J)*2*Z2/(ZW2) 

G  (  J  ;  =  G  (  J  )  *  2  *  Z  2  /  (  Z  1  *  Z  2 
F  <  .J  :  F  (  .)  )  *  Z  1  /  Z  2  ♦  W  i  J  i  /  •  W  t  1  «  C  <  1 

END  DO 

B  B  B  P  *  B  E  T  A  <  1 
AB -SMALLA»BETA  1 ! 

W  1C  T-  W  I  1  *C<  1  ' 

WPWC  =WP/W [C 1 
I  F  ■  SCS  . GT .  0  )  THEN 

CALL  SPACE'  AB , BB , SCS , R , EXP  1  , CB 


AVERAGED  IMPEDANCE’ 
’GAIN  PARAMETER' 
VELOCITY  PARAMETER' 
CIRCUIT  LOSS  (  <1 B  /  in 

LOSS  PARAMETER' 


X.  Si 


WPWC-0 
END  IF 
WQ  =  R  *W  P 

NORMZ=WlCl/UZERO 

YFINAL=ZF*NORMZ 

PRNT=DPC*NORMZ 

H-Z* NORMZ 

RETURN 

END 


$  EM A  // 

SUBROUTINE  ATTE N i  Y , F I RST .  I  AS  TOP ) 

C  ACCOMODATES  FOR  LOSS  DUE  TO  AN  ATTENUATOR,  AND  CALLS  THE 

C  DRIFT  TUBE  EQUATIONS  FOR  A  SEVER 

C  - 

IMPLICIT  DOUBLE  PRECISION  (A-H.OZ) 

INTEGER  SCS, SWITCH, SWLIST.SETU 

DOUBLE  PRECISION  I COS , I FUNCT , IS  I N ,  I  0 , K 1 , K2 , K3 , LA , LMAX , LOGlOE , 
&LOSS , MAG , NORMZ 

COMMON  A( 3 )  , ANG( 3 )  , AXIMP(3) , AXIMPN( 60, 3) , A0( 3 ) , B(3)  ,  BETA! 3 )  , 

&  B  E  T  A  0  (  3  ) ,  B  P  ,  C  (  3  ) ,  C  B  (  5  )  .CKTRADI60) ,0(3) ,DA(3) ,DPC,0PDB,DTHETA(3) , 
A0F(  3 )  ,  EBEAM, EBLNC , EM, EPS  I , ESAT, ETA , EWAVE , EXPl , F( 3 ) , FI LL, G( 3 ) , 

&G  AMA ( 3 )  , GS AT ,  I COS ( 3  5  ,  I FUNCT (64 ) ,  INTSKP,  ISIN(3) ,  I0.K1.K2.K3, 
fcLA( 20 , 6 ) , MAG( 3 ) , NATTN, NCHNG, ND , NE LEC T , NORMZ , NPR , NUMF , 0 LD PD B , 0 LDZ , 
fS.PDB,PDBl,PHI(64,3>,PHVEL(3),PHVELN(60,3),PI,PINI(3),PRNT,PSAT,PW, 
&RFLSS(3),RFLSSN(60,3),SCS,SETU,SKPINT, SMALL A , SMALLC,SWLIST(140), 

&  THE  T  A ( 3 )  , TPI , U  I 64 )  , UTERM ( 64)  , UZERO, VP( 3) , VREL, V0,W( 3)  , WP, WPWC , 

& Y FINAL, YH ALT, Y PRINT, Z , Z F , Z S AT , ZS TOP ( 1 40 ) 

LMAX  =  170  . 

LOGlOE- LOG  1  0 ( EXP(  I  .  )  1 

N  ATTN  =  N  ATTN  1 

DO  NSEC- 1,20  '20  ATTENUATOR  SECTIONS! 

ND-ND+  1 

YHALT -ZS TOP < ND )* NORMZ 

DO  J  =  1  , NUMF  'RECOMPUTE  LOSS  PARAMETER' 

LOSS =RFLSS(J)+SQRTfW(J)/W(  1  )  ) *  LA ( NSEC , NATTN ) 

D(  J  )  =LOSS*  .  05/ (  Wf  .!  ) /UZERO*C(  J  )  *LOG10E  )  'LOSS  PARAMETER' 

END  DO 

IF  ( LA ( NSEC , NATTN 1  .GT.  LMAX)  THEN  ! TREAT  AS  A  SEVER' 

DO  I- 1, NUMF 

A (  I  )  - M A X ( A 0 (  I  )  * . 0  1  DO , A (  I  ) * . 0 1  DO ) 

THE T A (  I  >  =  0 
G  (  I  )  -  0 
F(  I )  -  0 
END  DO 

CALL  DRIFT t  Y  ) 

Z  =  Y / NORMZ 

CALL  OUTPUT  f  F  I  RST 

IF  <  I.Ai  NSEC*  1  ,  NATTN  1  .  LE  .  LMAX)  THEN  'INITIAL  CONDITIONS’ 

DO  1=1, NUMF 
J  =  I 

CALL  FINTf  PHI  ,  U,  C  (  1  ;  ,  .1  ,  TS  I  N  ,  TCOS  ,  NKI.ECT  ,  UTERM  ,  P  I  ) 

ISINi  I  )  -TSIN 
ICOS (  II  TCOS 
END  DO 

W1C1  W  1  i  *  C  I  1  ! 

DO  N  I,  NUMF  'SET  INITIAL  CONI).  FOR  END  OF  SEVER' 

WNW  1  -  W  '  N  )  '  V*  (  1  ) 

WNW2  WNW 1 *WNW 1 


LMAX  ) 


! TREAT 


SEVER ' 


* . 0  IDO , A (  I  ) * . 0 1  DO ) 


LMAX  j 


INITIAL  CONDITIONS’ 


o  o 


GSAT=PDB1-PINI( 1 ) 

ESAT=ETA*100. 

ZSAT=Z 
END  I  F 

DTDZ  =  F( 1  )/NORMZ 

IF  ( ( NPR  .EO.  2)  .OR.  ( NPR  .EQ.  4)) 

&  CALL  ENERGY! A, AO ,C , U , NUMF , NELECT , EWAVE , EBEAM, EBLNC ) 

IF  ( FIRST  .NE.  I)  THEN 
WRITE! 98 , 5 ) 

IF  ! NPR  .EQ.  2  .OR.  NPR  .EQ.  4)  WRITE(98,10) 

IF  ( NPR  .EQ.  3  .OR.  NPR  .EQ.  4)  WRITE(98,20) 

5  FORMAT! /5X,  ’Z’  ,7X,  ’ PD  B ’  ,7X,  ’B’  ,5X,  ’ THETA ’,5X, ’DPDB’, 6X, ’DTDZ’, 
&  6X , ’ETA’ , 6X , ’ D_ ’ ) 

10  FORMAT! /IX ’EBE AM’ , 8X, ’EWAVE’, 8X, ’EBLNC’) 

20  FORMAT!/’  PHI(I,1)  ELECTRON  PHASE  ARRAY  (FUNDAMENTAL)’/ 

A  ’  2*C(1)*U(I)  NORM.  AC  VELOCITY  ARRAY’) 

F I RST= 1 
END  IF 

WRITE! 98, 25)Z, PDB , B ( 1 ) , THETA! 1) , DPDB , DTDZ , ETA  ,  D  (  1 ) 

25  FORMAT!/iX,F7.3,F9.3,F8.2,F9.3,F10.3,F10.4,F8.4,F9.4) 

IF  (NUMF  .GT.  1)  THEN 
DO  J- 2 , NUMF 

PW  =  2.  *C(  1 )  *I0*VREI,*A!  J  )  *A<  J  )  *(  1  .  -C!  J  )  *F(  J  ))  /(  1+C(  J  )  *B!  J  )  ) 

PD  B  =  1 0  . *  LOG  1 0 ( ABSIPW)  ) ^30  - 
PDBI(J)=PDB 

WRITE(98,30)PDB,B!J) , THETA!  J) 

30  FORMAT! 8X, F9. 3 , F8. 2, F9. 3 ) 

END  DO 
END  IF 

WRITE! 97, 33) Z, ! PDB I! N) ,  N=1,NUMF) . (THETA! Ml , M= 1 , NUMF ) 

WRITE! LU, 33 )Z , ! PDB I ! N) , N=1 , NUMF) 

IF  ! NPR  .EQ.  4)  WRITE(96,32)Z, (PHI! I, 1) , 1=1, NELECT) 

32  FORMAT! IX, F4. 3,64! IX, F5. 2) ) 

33  FORMAT! IX, F7 . 4 , 6!2X, F8. 4 ) ) 

IF  ((NPR  , EQ .  2)  .OR.  (NPR  .EQ.  4))  THEN 
WRITE (98, 35) EBE AM, EWAVE.EBLNC 
35  FORMAT! I  X , 3 ( G  1  2 . 5 , 2 X  )  ) 

END  IF 

IF  ((NPR  .EQ.  3)  .OR.  (NPR  .EQ.  4))  THEN 
WRITE (98,45 ) (PHI( I, 1) . 1=1, NELECT) 

WRITE (98, 45) ! 2. *C( 1 ) *U( I ) , 1= 1 , NELECT) 

45  FORMAT (4(1X,G12.4,3X)) 

END  IF 
RETURN 
END 

C  =  =  =  =  =  =  =  =  i  =  =  =  =  =  =  =  ==  =  ==  =  =  =  =  =  =  =  =  =  =  =  =  =  =  -  =  =  =-  =  =  =  =  r-  =  =  r.  —  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =;::: 

$EMA  // 

SUBROUTINE  OUTRS(Y) 

PRINTS  OUT  RESTART  DATA 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS , SWITCH, SWLIST , SETU 

DOUBLE  PRECISION  I  COS,  IFUNCT,  I  SIN,  10,  IOGB,  I1GB,K0GB,K1GB,K1.K2,K3, 
ALA, LT , MAG , NORMZ 

COMMON  A(3) , ANG(3) , AX  IMP ( 3 ) , AX IMPN! 60 , 3 ) , AO ! 3  > , B ( 3 ) , BETA ( 3 )  , 
ABETAO! 3) , BP,C( 3) ,CB(5 ) ,CKTRAD(60 ) , D( 3 ) , DA( 3 ) , DPC, DPDB, DTHETAi 3  , 
ADF<3) , EBE AM, EBLNC, EM, EPS  I , ESAT, ETA, EWAVE, EXP1 , F(3) , FILL, G( 3  )  . 

AG AM A! 3 ) , GS AT ,  I COS ( 3 ) , I FLNCT ! 64 ) ,  I  NTS KP , I S I N ( 3  )  ,  1 0 , K I , K2 , K3 , 
ALA(20,6) , MAG ( 3 ) , N ATTN , NCHNG , ND , NE LECT , NORMZ , NPR , NUMF , OLD PDB , OLD Z , 
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fcPDB , PDB 1 , PHI ( 64 , 3 ) , PH  VE  L I  3 ) ,PHVELN(60,3) ,PI,PINI(3) ,PRNT,PSAT,PW, 
4RFLSS(3) , RFLSSN(60 . 3  )  , SCS , SETU, SKPINT, SMALLA, SMALLC, SWLISTf 140 ) , 
&THETA ( 3 )  ,  TP  I , U ( 64 ) , UTERM ( 64 ) , UZERO , VP ( 3 ) , VRE L , VO , W ( 3 ) , WP, WPWC , 

& YF I N A L  , YHALT, YPRINT, Z , ZF, ZSAT, Z  S  TOP ( 140 ) 

WRITE ( 98 , 10) Y 

10  FORMAT ( / / ’  RESTART  INFORMATION  FOR  Y -  ’,F7.3/) 

WRITE(98,  20  ) 

20  FORMAT ( 6X ,  'A’  , 1 IX ,  ’DA’  ,8X,  ‘THETA’  ,7X,  '  DTHET A  '  ) 

DO  1=1, NUMF 

WRITE (98, 30) A(  I  )  . DA( I ) , THETA(  I )  , DTHET A  I  I ) 

30  FORMAT(4( IX, E10.4, IX)  ) 

END  DO 

WR I TE ( 98 , 40 ) 

40  FORMAT ( / ’  ELECTRON  PHASE  POSITIONS’) 

DO  I = 1 , NE LECT 

WRITE(98,30)(PHI(I,J),J=1, NUM  F I 
END  DO 

WRITE(98, 50) 

50  FORMAT ( / '  ELECTRON  VELOCITIES’) 

WRITE (98, 30) (U( I) . 1=1 , NELECT) 

IF  (SCS  .EQ.  2)  THEN 

WRITE (98, 60)  (CB( I )  ,  1  =  1 ,5) 

60  FORMAT ( / ’  COEFF.  FOR  POLYNOMIAL  EXPANSION  OF  SPACE  CHARGE’/lX, 
*.  5  (  E  1  0 . 4 , 2  X  )  ) 

END  IF 
RETURN 
END 


SUBROUTINE  E NE RG Y( A, AO, C.U, NUMF, NELECT, EWAVE, EBEAM, EBLNC) 
C  CHECKS  FOR  CONSERVATION  OF  ENERGY 

c  - 

IMPLICIT  DOUBLB  PRECISION  (A-H.O-Z) 

DIMENSION  A (NUMF) , AO (NUMF) ,C(NUMF) , U( NELECT) 

EMA  A, AO ,C, U, NUMF, NELECT, EWAVE , EBEAM, EBLNC 
EW A  VE  =  0 
EBEAM=0 
DO  1=1, NUMF 

EWA VE  =  A ( I ) * A ( I ) -AO ( I ) *  AO ( I ) +EWAVE 
END  DO 

EWA VE  =  2 . *C( 1 ) *  EW  A  V  E 
DO  1=1, NELECT 

CTERM=  1  . *2 .  *C  t  1  )  *U(  I ) 

EBEAM=CTERM*CTERM+EBE AM 
END  DO 

EBE AM  =  EB  E AM/NE  LECT -  1  . 

EBLNC =ABS ( EBEAM ) -ABS ( EWAVE ) 

EBLNC=EBLNC/( ABS( EBE AM )+ABS( EWAVE)) 

RETURN 

END 

C==-  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =-  =  =  -  =  r  =  “  =  =  ;  =  =  =  =  :rr-:r  “  :  :•  r 

SUBROUTINE  FI  NT (PHI, U, Cl, J.TSIN.TCOS, NELECT, UTERM, PI) 

C  CALCULATES  THE  BEAM  CHARGE  DENSITY 

r  _ _ _ _ _ _ _ _ _ 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DIMENSION  PHI ( 64 , 3  i  , U( NELECT) , UTERM f  NELECT 1 
EMA  PH  I, U, Cl, NELECT, UTERM. PI 
TWOP 1=2. *PI 
TS I N  =  0 
TCOS  =  0 


o  o 


DO  I = 1 , NE  LECT 

UTERM ( I ) = 1 . +2 . *C 1 *U( I ) 

TSIN  =  TSIN-*-SIN(PHI( I , J ) )/UTERM( I ) 

TCOS=TCOS+COS( PHI ( I , J ) ) /UTERM ( I ) 

END  DO 

RN  =  F  LO AT ( NELECT) 

TSIN=TSI N*TWOP I / RN 
TCOS - TCOS  * TWOP I / RN 
RETURN 
END 

C=========== rzr====r============r=====r=z====r-=-3-r-r:-r-- 

$EMA  /TEMP/ 

SUBROUTINE  INTI (PHI , UTERM, EXP  1 .NELECT,  IFUNCT, PI ) 
SOLVES  EXACT  SPACE  CHARGE  INTEGRAL  FOR  FORCE  EQUATION 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  I FUNCT ( NE LEC T ) ,  I  FUNSM  ,  PHI  (  64  ,  :i)  ,  UTERM  (  NE  LECT  ) 
EM A  PHI , UTERM, EXP  1 , NELECT,  IFUNCT, PI 
COMMON  /TEMP/  FTEMP(64,64) 

TWOP I  =  2  *P I 
NM1 -NELECT-1 
DO  1=1, NM 1 
IP1  =  1+  1 

DO  J  =  I P  1  , NELECT 

ZZ  =  PHI ( 1 , 1 ) -PH I  ( J ,  1  ) 

IF  ((ZZ  .GE.  TWOPI)  .OR.  (ZZ  .LT.  0)) 

8.  ZZ  =ZZ- TWOP  I  *  (  I  NT  (ZZ/ TWOP  I  )-(l.-SIGN(].D0,ZZ))*.5) 

IF  ( :( ZZ  .GE.  TWOPI)  .OR.  (ZZ  .LT.  0))  GO  TO  500 
S I NZZ  =  S I N ( ZZ  ) 

COSZZ=COS(ZZ  ) 

TERM 1=SINZZ/(EXP1-C0SZZ) 

TERM2  =  SINZZ/(  EXPUEXP1  -COSZZ  ) 

FTEMP ( I , J ) = ( ( P I -ZZ ) * . 5-2 . *ATAN ( TERMl )  +  AT AN ( TERM2 ) ) /TWOPI 
FTEMP ( J , I ) = -FTEMP ( I , J ) /UTERM ( I ) 

FTEMP  (  I  ,  J  )  =  FTEMP  (  I  ,  J  ) /'UTERM  (  J  ) 

END  DO 
END  DO 

DO  1=1, NELECT 
I FU  N  SM  =  0 
DO  J--1,  NELECT 

IFUNSM=IFUNSM+FTEMP( I , J  ) 

END  DO 

I FUNCT (  I)  =  I FUNSM*TWOP I / F  LO AT ( NM 1 ) 


END  DO 

500 

RETURN 

END 

C-  =  =  = 

$  EM  A 

/TEMP/ 

SUBROUTINE  INT2(PHI,UTERM,CB, NELECT, I FUNCT , P I) 

C  SOLVES  THE  POLYNOMIAL  APPROXIMATION  TO  THE  SPACE  CHARGE 

C  INTEGRAL  FOR  THE  ELECTRON  FORCE  EQUATION 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

DOUBLE  PRECISION  CB  (  5  )  ,  I  F  UNC  T  (  NE  LEC  T  )  ,  I  FUNSM,  PH  I  (  64  .  .'IS  , 
AUTERM ( NE i EC T  ) 

EM A  PHI , UTERM. CB , NELECT , I FUNCT , P I 
COMMON  /TEMP/  FTEMP ( 64 , 64  ) 

TWOPI -  2 . * P I 
NM1=NELECT  1 


DO  1=1, NM 1 
IP1=I+1 

DO  J= IP  1 , NE  LECT 

ZZ=PHI ( I, 1 ) -PHI ( J , 1) 

IF  ((ZZ  .GE.  TWOPI)  .OR.  (ZZ  .LT.  0)) 

.  ZZ=ZZ-TWOPI*( I NT (ZZ/ TWOPI) -( 1 . - S IGN ( 1 . DO , ZZ ) ) * . 5 ) 

IF  ((ZZ  .GE.  TWOPI)  .OR.  (ZZ  .LT.  0))  GO  TO  500 
ZZ=ZZ/PI-1 . 

FTEMP ( I , J ) =CB ( 3 ) *ZZ*ZZ*ZZ*ZZ*ZZ+CB ( 2 ) *ZZ*ZZ*ZZ+CB  (  1 ) *ZZ 
FTEMP ( J , I ) =-FTEMP( I ,J)/UTERM( I) 

FTEMP ( I ,J)=FTEMP(I,J)/UTERM(J) 

END  DO 
END  DO 

DO  1=1, NELECT 
I FUNSM=0 
DO  J=l, NELECT 

IFUNSM=IFUNSM+FTEMP( I, J) 

END  DO 

I FUNCT( I ) =IFUNSM*TWOPI/FLOAT(NMl ) 

END  DO 
RETURN 
ND 


$  EMA  // 

SUBROUTINE  DRFTEQ( Y , VALUES , EONS ) 

C  PHASE  AND  FORCE  EQUATIONS  FOR  THE  DRIFT  TUBE.  B ( N )  REPLACES 

C  DTHETA  IN  THE  PHASE  EQUATION  TO  INSURE  THE  CORRECT  INITIAL 

C  PHASE  FOR  THE  CIRCUIT  WAVE  IN  THE  NEXT  SECTION 

C  DPHI( I , N ) — PHASE  EQUATION;  DU(I)--FORCE  EQUATION 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS , SWITCH, SWLIST, SETU 

DOUBLE  PRECISION  ICOS , I FUNCT , I S I N , 1 0 , K 1 , K2 , K3 , LA , MAG , NORMZ , 
ftDU(64)  ,  DPHK64.3)  ,  EONS  (  268)  ,  VALUES  (  268  ) 

COMMON  A ( 3 ) , ANG ( 3 ) , AXIMP(3) , AX IMPN ( 60 , 3 ) , A0(3) , B(3) , BETA ( 3 ) , 
&BETAO ( 3 ) , BP , C ( 3 ) , CB ( 5 ) ,CKTRAD(60) , D(3 ) , DA(3) , DPC , DPDB , DTHETA ( 3 ) , 
&D  F ( 3  ) , EBEAM, EBLNC , EM, EPS  I , ESAT, ETA, EWAVE, EXP1 , F(3),FILL,G(3), 

&G  AMA ( 3 )  , GS AT ,  ICOS(3)  , I FUNCT (64) , INTSKP, IS  IN ( 3 ) , 10 , K1 , K2 , K3 , 

&LA ( 20 , 6 ) , MAG ( 3 ) , NATTN , NCHNG , ND , NELECT , NORMZ , NPR , NUMF , OLDPDB , OLDZ 
&PDB , PDB 1 , PH 1(64, 3) , PHVEL(3) , PHVELN(60 ,3),PI,PINI(3), PRNT, PS AT, PW 
&RFLSS ( 3 ) , RFLSSN ( 60 , 3 ) , SCS , SETU , SKP I  NT , SMA LLA , SMALLC , SWLI ST ( 140) , 
&THETA( 3) , TP  I , U(64) , UTERM(64) , UZERO, VP( 3 ) , VREL, VO, W(3) , KP, WPWC , 
&YFINAL, YHALT, YPRINT, Z ,ZF,ZSAT,ZSTOP(140) 

DO  1=1, NELECT  'INITIAL  VALUES  FOR  DRIFT  EQNS ! 

U (  I  )  =VALUES( I ) 

END  DO 
DO  J= 1 , NUMF 

DO  1=1, NELECT 

PHI(I,J)=VALUES( 1+ J* NELECT ) 

END  DO 
END  DO 

IF  ( S  KP  T  NT  .EQ.  0)  THEN  'COMPUTE  SPACE  CHARGE  FIELDS  ft  INTEGRAL 
DO  1=1, NUMF 

J=  I 

CALL  FINTIPHI, U,C( 1 ),J,TSIN, TCOS , NELECT, UTERM, PI ) 


INTEGRA  L 


CALL 
END  DO 
IF  ISCS 
IF  ISCS 

END  IF 


NELECT, UTERM, PI ) 


CALL  I NT2 ( PHI , UTERM, CB , NELECT , I FUNCT, PI) 
CALL  INTI  fPHI,  UTERM, EXP1  .NELECT,  IFUNCT.PI  ) 


o  n 


SKPINT=SKPINT+1 

IF  ( SKPINT  .GT.  INTSKP )  SKPINT=0 

DO  1=1 , NELBCT  ! DRIFT  TUBE  EQUATIONS' 

CTERM= 1 . +2. *C( 1 ) *U(  I  ) 

DO  N= 1 , NUMF 

BTERM= 1 .  +  C  (  N  )  *  B  (  N  ) 

DPHI( I,N)=B(N)+2. *W(N) *U( I )/(W( 1) #CTERM ) 

DU ( I ) =WPWC*WPWC  *  I FUNCT (  I) / ( BTERM*CTERM) 

END  DO 
END  DO 

DO  1=1, NELECT  !STORE  DRIFT  EQUATIONS  IN  ARRAY  EQNS ! 

EQNS ( I )  =DU(  I  ) 

END  DO 
DO  J  =  1 , NUMF 

DO  1=1, NELECT 

EQNS ( I+J*NELECT) =  DPH  I  (  I ,  J  1 
END  DO 
END  DO 
RETURN 
END 

$EMA  // 

SUBROUTINE  DRIFT(Y) 

C  SETS  THE  INITIAL  CONDITIONS  AND  CALLS  THE  DIFFERENTIAL 

EQUATION  SOLVER  FOR  THE  DRIFT  TUBE  EQUATIONS 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS , SWITCH, SWLIST , SETU 

DOUBLE  PRECISION  ICOS ,  I FUNCT , IS  IN , 10 , K1 , K2 , K3 , LA , MAG , NORMZ , 
&VALUES( 268) 

COMMON  A ( 3 ) , ANG ( 3 )  , AX  IMP ( 3 ) , AX IMPN ( 60 , 3 ) , AO (3) ,B(3) , BETA  I  3 )  . 
&BETA0I3)  , B  P , C ( 3 ) ,CB(5)  ,CKTRAD(60) , D(3) ,DA(3)  , DPC , DPDB , DTHETAI 3  )  , 

& D F ( 3 ) , EBEAM, EBLNC, EM, EPS  I , ES AT , ETA , EWAVE , EXP1 , F ( 3 ) , F I LL , G ( 3 ) , 

&G  AMA ( 3  )  , GS AT , ICOS ( 3 ) , I FUNCT (64  ) ,  INTSKP,  I S I N ( 3 )  ,  I  0 , K 1 , K2 , K3 , 

ALA (20, 6) , M AG ( 3 ) , NATTN , NCHNG , ND , NELECT , NORMZ , NPR , NUMF , OLDPDB , OLDZ , 
&PDB.PDB1  , PHI ( 64 , 3) ,PHVEL(3) , PH VE LN ( 60 , 3 ) ,PI,PINI(3) , PRNT , PS AT , PW , 
&RFLSS(3),RFLSSN(60,3)  , SCS , SETU , S KP I  NT , SMALL A , SMA LLC , SWL I S T i  140)  , 
&THETA ( 3 ) , TP  I , U(64 ) , UTERM(64) , UZERO, VP(3) , VREL, VO, W( 3 ) , WP, WPWC , 
&YFINAL, YHALT, YPRINT, Z , Z F , ZS AT , ZSTOP (  1 40  ) 

EXTERNAL  DRFTEQ 

DO  1=1, NELECT  IINITIAL  VALUES  FOR  DRIFT  EQNS' 

VALUES (  I ) =  U ( I ) 

END  DO 
DO  J  =  1 , NUMF 

DO  1=1, NELECT 

VALUES! I+J*NELECT)=PHI( I, J) 

END  DO 
END  DO 

NEQN=NUMF*NELECT+NELECT 
ABSERRr 1 . D-09 
I  FLAG  =  1 
RE  LE  RR  =  1  . D - 0  3 

CALL  DESOLVf  DRFTEQ , NEON, VALUES, Y,YH ALT, RE LERR.ABS ERR,  IFLAC, SETU) 
DO  1=1, NELECT 

U( I ) = VALUES!  I  ) 

END  DO 
DO  J  =  1 , NUMF 

DO  1=1, NELECT 

PHI ( I ,J)=VALUES( I+J*NELECT) 


END  DO 
END  DO 
RETURN 
END 

$ EMA  // 

SUBROUTINE  DIFFEQ(Y, VALUES, EQNS) 

C  ROWE’S  LARGE  SIGNAL  EQUATIONS 

C  D  A  (  N  )  8.  DG  I  N  )  :  CIRCUIT  WAVE  AMPLITUDE  EQUATION 

C  DTHETA(N)  &  D  F ( N ) :  CIRCUIT  WAVE  PHASE  EQUATION 

C  L?HI(I,N):  PHASE  RELATION  EQUATION 

C  D U ( I )  :  ELECTRON  FORCE  EQUATION 

C  - 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS , SWITCH , SWLIST, SETU 

DOUBLE  PRECISION  ICOS , I FUNCT ,  I S I N ,  I  0 , K 1 , K2 , K3 , LA , MAG , NORMZ , 

&DU ( 64 ) , DPHI ( 64, 3) , DG ( 3 ) , EQNS ( 268 ) , VALUES(268) 

COMMON  A ( 3 )  , ANG(3  )  , AXIMP(3) , AXIMPNI60, 3)  , AO ( 3  )  , B (3 ) , BETA( 3  )  , 
&BETA0 ( 3 ) , BP , C ( 3  ) ,CB(5)  , CKTRAD (60 ),D(3),DA(3),DPC,DPDB, DTHET A ( 3 ) , 
&DF(3) , EBEAM, EBLNC, EM, EPS  I , ESAT, ETA, EWAVE, EXP1 , F ( 3 ) , FILL,G( 3) , 

&G AMA ( 3 )  ,GSAT,  ICOS (3) ,  I  FUNCT (64) ,  INTSKP, ISIN(3),I0,K1,K2,K3, 

&LA ( 20 , 6 ) , MAG ( 3 ) , N ATTN , NC HNG , ND , NE LECT , NORMZ , NPR , NUMF , 0 LDPD B , OLDZ , 
&PDB , PDB1 , PHI (64 , 3) , PHVEL(3) , PHVELN(60, 3) , PI , PINI(3) , PRNT, PS  AT, PW, 
&RFLSS(3) , RFLSSN(60, 3)  , SCS , SETU , S KP I  NT , SMALL A , SMAL LC , SWL I S T ( 140) , 
&THETA( 3 ) , TP  I , U( 64 ) , UTERM(64) , UZERO, VP( 3 ) , VREL, VO , W( 3 ) , WP , WPWC , 
&YFINAL, YHALT, YPRINT, Z , ZF, ZSAT, ZST0P( 140) 

DO  1  =  1 , NELECT  (INITIAL  VALUES  FOR  DIFF.  EQNS! 

U( I)= VALUES! I) 

END  DO 
DO  J-l.NUMF 

DO  1=1, NELECT 

PHI(  I,  J)  =  VALUES  (  I  +- J*NELBCT  ) 

END  DO 
END  DO 

K  =  NE  LEC  T* I NUMF  + 1 ) 

DO  N=  1  , NUMF 
J  =  N  +  K 

G ( N ) =VALUES(J) 

F ( N ) =VALUES( J+NUMF) 

A(N) =  V  A  LU  E  S ( J+2*NUMF ) 

THETA  (N)  =  VALUES  (J-t- 3* NUMF) 

END  DO 

IF  ( SKPINT  .EQ.  0)  THEN  !  COMPUTE  SPACE  CHARGE  FIELDS  «.  INTEGRAL! 
DO  1=1, NUMF 
J=I 

CALL  FINTf  PHI , U,C( 1 ) , J, TSIN, TCOS , NELECT, UTERM, PI) 

I S  I  N ( I ) =TSIN 
ICOS (I ) =TCOS 
END  DO 

IF  (SCS  .EQ.  2)  CALL  I NT2 ( PH  I , UTERM , CB , NE LECT ,  I FUNCT , P I ) 

IF  (SCS  .EQ.  1)  CALL  I NT  1 ( PH I , UTE RM , E XP 1 , N E LECT ,  I  FUNCT , P  I  ) 

END  IF 

SKPINT=SKPINT*I 

IF  (SKPINT  . G  T .  INTSKP)  SKPINT  =  0 
W1C1 =W( 1  )  *C(  1) 

DO  N= 1 , NUMF 

WNW1 =W( N I /W( 1 ) 

WNW2  =  WNW 1 *WNW1 

BTERM=  (  l.  +-C  (  N  )  *  B  (  N  )  )/C(  1  ) 


DTERM=2 . *C(N)*D(N) 

FTERM=F(N)-W(N)/W1C1 

ZN0RM=AXIMP(N)/AXIMP(1) 

D A ( N ) =G ( N )  ! D IFFERENTIAL  EQUATIONS! 

DTHETA( N  )  =  F  (  N  ) 

DG(N) =  A(N) *( FTERM**2-WNW2*BTERM**2)-ZNORM*WNW2*BT£RM/PI* 

&  ( ICOS( N) +DTERM*ISIN( N) ) 

DF(N)=DTERM*WNW2*BTERM**2-(2.*G(N  )  *  FTE  RM  +  Z  NORM*WNW2  *  BTERM/ P  I  * 
&  (  ISIN(N)-DTERM*ICOS(N) ) )/A(N) 

DO  I=1,NELECT 

CTERM= l . +2 . *C (  1 )  * U (  I  ) 

DPHI ( I , N) =2. *WNWI*U( I ) / CTERM- FIN) 

END  DO 
END  DO 

DO  I = 1 , NE LECT 

CTERM= 1  . +2 . *C ( 1 ) *U(  I  ) 

D  U ( I ) =0 
DO  N=  1  , NUMF 

FTERM=  F ( N ) -W(N)/W1CI 
BTERM=1 . +C( N) *B( N) 

CS=COS( PHI ( I , N) ) 

SN  =  S  INI  PHI ( I , N)  > 

DU( I )  =  DU  ( I )+(G(N)*CS+A(N)*FTERM*SN) *C( 1 ) 

DU( I )  =  (DU( I ) +WPWC  *WPWC* I FUNCT ( I ) /BTERM) /CTERM 
END  DO 
END  DO 

DO  1=1, NELECT  'STORE  DRIFT  EQUATIONS  IN  ARRAY  EONS' 

EQNS (  I  ) =DU( I ) 

END  DO 

DO  J  -  I  , NUMF 

DO  1=1, NELECT 

EQNS ( I+J*NELECT) -DPHI ( I, J) 

END  DO 
END  DO 

K  =  NELECT*(NUMF»-1  ) 

DO  N= 1 , NUMF 
J  =  N  +  K 

EQNS(J)=DG(N) 

EQNS(J+NUMF)=DF(N) 

EQNS( J+2*NUMF) =DA(N) 

EQNS(J+3*NUMF)=DTHETA(N) 

END  DO 
RETURN 
END 

C====================================:==r;=====-===^=======r==r====-== 

$  EMA  // 

SUBROUTINE  DIFF(Y) 

C  'SETS  THE  INITIAL  CONDITIONS  AND  CALLS  THE  DIFFERENTIAL 
C  EQUATION  SOLVER  FOR  ROWE’S  LARGE  SIGNAL  EQUATIONS 

C  - 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

INTEGER  SCS, SWITCH, SWLIST, SETU 

DOUBLE  PRECISION  ICOS ,  I  FUNCT ,  l  S  I  N , I  0 , K 1 , K2 , K3 , L A , M AG , NORMZ , 

& VA  LUES ( 268 ) 

COMMON  A ( 3 ) , ANG ( 3 ) , AX  IMP ( 3 ) , AX IMPN ( 60 ,  3  )  ,  AO  (  3  )  ,  B  (  3  )  , BETA ( 3 )  , 
ABETAO ( 3 ) ,BP,C(3) ,CB(5) , CKTRAD ( 60 ),D(3),DA(3),DPC,DPDB, DTHETA ( 3 ) , 
ADF(3) , EBEAM, EBLNC , EM, EPS  I , ESAT, ETA , EWAVE ,EXP1,F(3),FILL,G(3), 
AGAMA(3>  ,GSAT,  ICOS  (  3  )  ,  I  FUNCT  (64)  ,  INTSKP,  I S  I  N  (  3  )  ,  10, K1  ,K2,K3, 

ALA  (  20 , 6  )  ,  MAG  (  3  )  ,  NATTN,  NCIING  ,  ND,  NELECT,  NORMZ  , NPR , NUMF , OLDPDB . OLDZ 


o  o  o  o  o  o 


S.PDB  ,  PDR1  ,  PHI  (  64 , 3  ,  PHVKL!  3  )  ,  PHVELN!  60, 3)  ,  PI  ,  PINI  (3)  ,  PRNT,  PS  AT,  PW, 
&RFI.SS  (  3  )  ,  RFLSSN(  (30 , 3  )  ,  SCS  ,  SETU  ,  SKPINT,  SMALLA  ,SMALLC,SWLIST(140), 

S. THETA  (  3  )  ,  TP  I  ,  U(  6  4  >\  UTERMI  64  )  ,  UZERO  ,  VP(  3  )  ,  VREL  ,  VO  ,  W(  3  )  ,  WP  ,  WPWC  , 
&YFINAL, YHALT, YPRINT, Z , ZF, ZSAT, /STOP ( 140) 

EXTERNAL  DIFFEQ 

DO  I-1,N ELECT  'INITIAL  VALUES  FOR  DIFF  EONS' 

VALUES! I } - U ( I > 

END  DO 
DO  J  =  1  , NUMF 

DO  1=1, NELECT 

VALUES!  I-*-J*NELECT)  -PHI!  I  ,  J  ) 

END  DO 
END  DO 

K-NELECT*!NUMF+1 ) 

DO  N  = 1 , NUMF 
J  =  N  +  K 

VALUES(J)=G!N) 

VALUES! J+NUMF)=F(N) 

VALUES(J+2*NUMF)-A(N) 

VALUES ( J+3*NUMF) = THETA! N) 

END  DO 

NEQN=4*NUMF+NELECT*NUMF+ NELECT 
AB  S  E  R  R  r 1 . D-09 
I  F  L AG  =  1 
RELERR: I . D-03 

CALL  DESOLV! D I FFEQ, NEON, VALUES , Y, YPRINT, RELERR, ABSERR,  I  FLAG , SETU) 
DO  1=1, NELECT 

U!  I  )  = VALUES ( I ) 

END  DO 
DO  J -  1  , NUMF 

DO  1=1, NELECT 

PH  I (  I  , J ) =  V  A  LUES ( 1  + J*NE  LECT ) 

END  DO 
END  DO 

K-NELECT*!  NUMF*1  ) 

DO  N  = 1 , NUMF 
J  -  N  +  K 

G(N)=VALUES(J) 

F( N) = VALUES ( J+NUMF) 

A! N)  = VALUES ( J+2*NUMF) 

THETA(N)=VALUES(J+3*NUMF) 

END  DO 
RETURN 
END 

C=-=========r=z========^====================-===::===========r======= 

$ EMA  /ABC/ 

SUBROUTINE  DESOL V ( F , NEON , Y, T, TOUT, RELERR, ABSERR, I  FLAG, SETU) 

C  INTEGRATES  UP  TO  268  DIFFERENTIAL  EQUATIONS  OVER  THE  INTERVAL 

FROM  T  TO  TOUT.  DESOLV  CALLS  THE  INTEGRATOR  STEP,  AND  THE 
INTERPOLATION  ROUTINE  INTRP.  ALL  THREE  SUBROUTINES  ARE 
LISTED  AND  FULLY  DOCUMENTED  IN  THE  TEXT  "COMPUTER  SOLUTION 
OF  ORDINARY  DIFFERENTIAL  EQUATIONS:  THE  INITIAL  VALUE  PROBLEM” 

BY  L.F.  SHAMPINE  &  M.K.  GORDON. 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

LOGICAL  START, CRASH, STIFF 
INTEGER  SETU 

DIMENSION  Y ! NEQN )  , PS  I ( 1 2 ) , Y Y ( 268 ) , WT ( 268 ) , PHI ( 268 , 1 6 ) , 
iP( 268 ) , YP ( 268 ) , YPOUT!268) 


k  n  fc.fl  Jil  AVkkV  J"k  .V  -A A/, 


-  61 


COMMON  /CDE/  YY , P , YP , PS  I , U , FOURU 

COMMON  /ABC/  PHI.WT.YPOUT 

EMA  SETU, TOUT 

EXTERNAL  F 

DATA  MAXNUM/2000/ 

EPS  =  MAX ( RELERR, ABSKRR ) 

I S  N  =  ISIGN( 1 ,  I  FLAG ) 

I  FLAG- I ABS (  IFLAG) 

IF  ( SETU  .EO.  0)  THEN 
CALL  MACHIN(U) 

FOURU=4  .  *U 
SETU- 1 
END  IF 

IF  (NEON  .LT.  1  .OR.  NEON  .GT.  268)  CALL  BADIN 
IF  (T  .EQ.  TOUT)  CALL  BADIN 

IF  (RELERR  .LT.  0.0  .OR.  ABSERR  .LT.  0.0)  CALL  BADIN 

IF  (EPS  .LE.  0  -OR.  IFLAG  .EQ.  0  )  CALL  BADIN 

IF  (IFLAG  .NE.  1  .AND.  T  .NE.  TOLD)  CALL  BADIN 

IF  (IFLAG  .LT.  1  .OR.  IFLAG  .GT.  5)  CALL  BADIN 

DEL-TOUT-T 

ABS DE  L  =  ABS ( DEL) 

TEND=T+ 10 . 0*D£L 
IF  (ISN  .LT.  0)  TEND -TOUT 
NOS  T  E  P  =  0 
K  LE4  =  0 

STIFF-.  FALSE. 

RELEPS-RELERR/EPS 
ABSEPS= ABSERR/EPS 

IF  (IFLAG  .EQ.  1  .OR.  ISNOLD  .LT.  0)  THEN 
STARTs. TRUE. 

X  =  T 

DO  L= 1 , NEON 
YY( L) =Y( L) 

END  DO 

DELSGN=SIGN( 1 . DO , DEL) 

H=SIGN(MAX( ABS ( TOUT- X) , FOURU* ABS ( X) ) , TOUT-X) 

END  IF 

10  IF  (  DELSGN*DEL  .  GT.  0.0  . AND .  ABS(X-T)  .GE.  ABSDEL )  THEN 
CALL  I NTRP ( X , YY , TOUT , Y , YPOUT , NEON , HOLD , PHI , PS  I ) 
IFLAG=2 
T  =  TOUT 
TOLD=T 
ISNOLD -  ISN 
RETURN 
END  IF 

IF  (ISN  .LE.  0  .OR.  ABS ( TOUT-X )  .LT.  FOURU* ABS ( X  )  )  THEN 
H^TOUT-X 
CALL  F ( X , Y Y , YP ) 

DO  L  = 1 , NEON 

Y(L)=YY(L)+H*YP(L) 

END  DO 
IFLAG  -  2 
T=TOUT 
TOLD=T 
ISNOLD^  ISN 
RETURN 
END  IF 

IF  (NOSTEP  .GE.  MAXNUM )  THEN 
I  FLAG  = I SN*4 
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SUBROUTINE  BADIN 

C  WARNS  USER  OF  BAD  INPUT  TO  DESOLV  AND  STOPS  PROGRAM 

EXECUTIO 

WRITE! LU, 10) 

10  FORMAT!’  INPUT  TO  THE  DIFFERENTIAL  EQUATION  SOLVER  1 

tS  BAD ' ’ ) 

STOP 

END 

SUBROUTINE  MACHIN(U) 

CALCULATES  THE  SMALLEST  POSITIVE  NUMBER  U,  SUCH  THAT 
1.0+U  >  1.0.  U  IS  COMPUTED  APPROXIMATELY  AS  A  POWER 
OF  1/2 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

HA  LFU  r . 5 

50  TEMPI  =  I  .  -t-HALFU 

IF  !  TEMP 1  . LE.  1 .  )  GO  TO  100 
HALFU= . 5*HALFU 
GO  TO  50 
100  U=2 . *HALFU 
RETURN 
END 

SUBROUTINE  S TE P(XtY,  F, NEON, H, EPS, WT, START, HO LD,K,KOLD, CRASH 


PHI.P, YP.PSI.U) 

INTEGRATES  THB  DIFFERENTIAL  EQUATIONS  OVER  A  STEP  (X  TO  X  +  H)  USING 
THE  MODIFIED  DIVIDED  DIFFERENCE  FORM  OF  THE  ADAMS  PECE  FORMLAS. 


IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

LOGICAL  START, CRASH, PHASE1 , NORND 

DIMENSION  Y ( NEON ) , WT ( NEQN )  ,PHI(NEQN,  16)  ,P(NEQN)  , YPfNEQN)  , 
I.PSI  (  12)  ,  GSTR  (  13  )  ,  TWO  (  1  3  )  ,  ALPHA  (  1  2  )  ,  BETA  (  1  2  )  .  S  IG  (  1  3  )  ,  W  (  l  2  >  , 
8>V(  12) ,G( 13) 

COMMON  /CSTEP/  GSTR , TWO , ALPHA , BETA , S IG , W , V , G 
EMA  PHI.WT 
EXTERNAL  F 

TWO( 1 ) =2. 

TWO ( 2 ) =4. 

TWO ( 3 ) =8 . 

TWO ( 4 ) = 1 6 . 

TWO( 5) =32 . 

TWO(6) =64. 

TWO ( 7 )  =  1 28 . 

TWO ( 8 ) =  256 . 

TWO ( 9 ) =5  12 . 

TWO ( 1 0 ) = 1 024 . 

TWO (  1 1 ) =2048 . 

TWO ( 12) =4096. 

TWO ( 1 3 ) =81 92  . 

GSTR ( 1 ) = . 5 
GSTR ( 2 ) = . 0833 
GSTR ( 3 ) = . 0417 
GSTR ( 4 ) = . 0264 
GSTR ( 5 ) = . 0 1 88 
GSTR ( 6 )  =  . 0 1 4  3 
GSTR(7)=. 0114 
GSTR(8)=. 00936 
GSTR  f  9 )  =  . 00789 
GSTR ( 10) = . 00679 
GSTR(11)=. 00592 
GSTR ( 12) =. 00524 
GSTR ( 13)=. 00468 
G  (  1  )  =  1 . 

G  (  2  )  =  .  5 
S  IG (  1)  =  1. 

CRASH= . TRUE . 

TWOU  =  2 . *U 
F  0  U  R  U  =  4 . *U 

IF  (ABS(H)  .LT.  FOURU*ABS ( X ) )  THEN 
H=SIGN(FOURU*ABS(X) ,H) 

RETURN 
END  IF 

P5EPS  = . 5  *EPS 
ROUND  =  0 . 

DO  L= 1 , NEQN 

ROUND  =  ROUND+(Y( L ) /WTf  L ) ) **2 

END  DO 

ROUND  =  TWOU*SQRTf  ROUND ) 

IF  (P5EPS  .LT.  ROUND)  THEN 
EPS=2. *ROUND*( 1 . +FOURU) 

RETURN 
END  IF 

CRASH= . FALSE . 


GO  TO  10 


IF  f  .  NOT .  START  > 

CALL  F  i  X  ,  Y  .  V  P  > 

S  UM  :  0 . 

DO  L -  1 , NEON 

PHI(  L,  1)  YP  L 
PHI( L, 2  0. 

SUM  SUM*  (  YP(  I.  WT(  I.  A  A  J 
END  DO 

S  UM  =  SORT f  SUM  , 

ABSH- ABS ( H 

IF  (EPS  .LT.  18. ASUMAHAH '  A  B  SH  .  2  5  *  SORT  K  PS  /  S  CM  ) 

H-SIGNi  MAX(  ABSH,  FOUR!)  *  ABS  i  X  ,H 
HOLD -0 . 

K  =  1 

KOLD  0 

START- . f ALSE . 

PHASE  1  -  . TRUK . 

NORND = . TRUE . 

IF  (P5EPS  .GT.  100. AROUND >  GO  TO  10 
NORND  -  .  FALSE . 

DO  L- 1 , NEON 

PHI  L,  15  >  -  0. 

END  DO 
I  F  A  I  L  -  0 
KP  1  -  K+  1 
KP2  =  K  *2 
KM  1  - K  -  1 
KM 2  - K  -2 

IF  i  H  ■  NE  .  HOLD  '  NS  0 
NS -M  I  NO •  NS* 1 , KOLD  *  1 
NSP1-NS+1 

IF  IK  . LT.  NS  1  GO  TO  50 
BETA ( NS ) = 1 . 

REALNS=NS 

ALPHA! NS ) -  1  . ,  REALNS 
TEMPI -H»REALNS 
SIGi NS PI) = 1 . 

IF  ( K  . GE .  NSP1 )  THEN 
DO  I  NS  PI, K 
I M 1  =  I  —  1 

TEMP2  PS  I !  I M 1  ) 

PSI  (  IM1  )  -  TEMPI 

BETA!  I  ;  -  BETA!  IMl  ,  APS  I (  IM1  :  TEMP 2 
TEMPI =TEMP2*H 
ALPHA I  I  >  -  H / TEMP  1 
R  E  A  L  I  =  I 

SIGI  I  ♦  1  )-REALI*ALPHA(  I  )  A  S  I G  <  I) 

END  DO 
END  IF 

PS  I  I K ) -  TEMP  1 
IF  INS  . LE .  1 )  THEN 

DO  I Q  z 1 , K 

TEMP3- 10* (  IQ*  1) 

VI  I  0 )  =  1 . / TEMP3 
W( IQ  )  -  V  I IQ) 

END  DO 
GO  TO  40 
BND  IF 

IF  I K  . LE .  KOLD )  GO  TO  30 


1 


TEMP4=K*KP1 
V  K  1 .  TEMP 4 
NSM2 -  NS  -  2 

I  p  nsm2  .  lt  .  i  <;o  ro  ;io 

DO  J  -  1  , NSM2 
I  K  J 

V  ;  I  =  V  I  A  ;  PH  A  I  »  1  *V ,  [  ♦  1 
END  DO 

LIMIT1  KP1  NS 
T  E  M  P  5  ALPHA  N" 

DO  10  l  ,  L I M !  r  1 

viiq  -v.ro  tempo  *  v  i  o  *  i  ■ 
w  i  r  o  v  i  o 
END  DO 
1 7  N  S  P  1  W  1 
NSP2  N  S  ♦  2 

IE  KP1  .OK.  NM'2  I  MEN 
DO  I  - N S  P 2 ,  K  P  1 
L  I  M  I  T  2  K  i '  2  I 
TF.MPH  All’  H  A  I  1 
DO  10  1  ,  I  I  M  I  I  2 

w  :  ; g  W  lo  r EMPH  *W  I  O  •  I 
END  DO 
0  I  w  ; 

END  DO 
END  I  E 

II-  r.  .  OK  N  .  I  HE  N 
lm  I  N..P  1  .  k 

TEMPI  BETA  r 
DO  !.  i  .  NEON 

PHI  I  ,  I  TEMPI  *  PH  I  1,1 
END  DO 
END  DO 
END  IE 
DO  1.  1  ,  NEON 

PHI  I.  ,  KI  2  PHI  I.  ,  KP  1  i 
PH  I  I.  ,  hP  i  ')  . 

P<  L '  0  . 

END  DO 
DO  I  1  ,  h 
I  k  P  1  I 
I  P  1  I  •  1 
TE  MP2  Oil: 

DO  I.  1  ,  NEON 

P  L  P  I  •  TEMP2  *PH  I  (  L  ,  I  ) 

PHI  I,  .  !  PHI  I.  ,  I  >  ♦  PH  I  (  L  ,  I  P  I 
END  DO 
END  DO 

IE  '  NOT.  NOHND  THEN 
DO  L  I, NEON 

T  A  I!  H  *  P  i  I.  ,  PHI  l  I.,  15  ) 

P '  I,  '  V  I  -TAP 
PH  I  (  I.  ,  1  0  i  P.  I  i  -  Y(  I.  )  )  T A II 
END  DO 
(70  TO  HO 
END  IE 
DO  L  1 , NEON 


P  II.  )  -  Y  (  I.  >  *  H  *  P  i  I.  ; 

END  DO 


X  O  L  D  X 
X  X  »  H 

ABSH  ABS  H) 

CAM.  F(  X  ,  P  ,  YP  ) 

EHKM2  0 . 

E  H  K  M  1  0  . 

E  W  K  0  . 

DO  I.  l.NEQN 

T  E  M  P  3  -  1 .  /WT (  L ) 

TEMP4^YP1 L ) -PHI ( l , 1 ) 

IP  ( KM2  )  90,80,70 

E  HKM2  =  ERKM2  + (  (PHI(L,KMl)*TEMP4)*TEMP3)*»2 
E  R  KM  1  -  E  R  KM  1  ♦•  (  (PHI(  L  ,  K  )  *T£MP4  )  * TEMP 3  )  **2 
E  K  K  ;  E  R  K  ♦ ( TEMP4*TEMP3 ) *  *  2 

E  N  D  DO 

I P  KM2 )  120,110,  100 

t  HKM2  ABSH*S IG( KM1 ) *GS  TR ( KM2 ) *S0RT( ERKM2 ) 

PR KM  1  -ABSH*SIG(K)*GSTH(KM1  )  *SQRT( ERKM1 ) 
TEMPS  ABSH*SQRT(  ERK  ) 

K  HR  r E M P 5  * .  G (  K )  G (  KP  1  )  ) 

FHK  T  F  M  P  fi  *  S  I  G  (  KP1  )  *GSTR(  K  ) 

KNEW  K 

I  ^  KM2  1  GO  ,  1  40 ,  1 30 

IP  MAX  ERKM 1 , ERKM2 )  .  LE .  E  HK  >  KNEW  KM  1 

go  ro  iso 

II  KHKMl  .LE.  . 5  *ERK )  KNEW  -  KM  1 

IP  ERR  LE .  EPS )  GO  TO  1 90 
PHASE  1  .  PA  l.SE  . 

x  x  u  i .  n 

I'O  i  i  ,  K 

TEMPI  1  .  / B  E  T  A (  I  l 
I  PI  I  ♦  1 
DO  L  1 , NEON 

PHI  L  ,  I  i  TEMP  1  *  (  PHI  (  L  ,  I  )  PHKL.IP1M 
E  N  f)  D  0 
END  D  0 

IP  K  . G  E .  2)  THEN 

DO  I  2 , K 

PS  I  I  1  >  -  PS  I f  I  )  -H 
END  DO 

END  IP 

I  P  A  I  L  I  P  A  I  L  ♦  1 
T  E  M  P  2  S 

IP  ‘  I  PA  I L  3)  180,170,160 

IP  '■  PS  EPS  .  LT .  .  2  S  *  E  RK  )  TEMP2 - SORT ( PS E PS / E 

KNEW  1 

H  -  J  E  M  P  2  *  H 
K  KNEW 

IP  'ABS(H)  LT.  FOURU*ABSfXM  THEN 
CRASH  .  TRI/E  . 

H  SIGN!  FO  U  R  U  *  A  B  S ( X  )  , H ) 

EPS  EPS  *  EPS 
RETURN 
END  IF 
GO  TO  20 
KOLD  -  K 
HOLD  H 

TEMP  1  H*G I KP1  i 
I P  ( NORND  )  GO  TO  200 


1  >  PS  I  E  I  ) 


180, 170, 160 

.  2  S  *  E  RK  )  TEMP2-SORT( PSEPS/ERK  ) 


FOURU*ABS ( X  I 


V*“- 


DO  L ' 1 , NEQN 

RHO= TEMPI  *(  YP(  L)  -PHI  (  L.  1)  )  -  PHI  (  L,  16 
Y ( L) =P( L ) ♦ RHO 
PHI ( L, 15 ) = ( Y( L) -P( L) ) 'RHO 
END  DO 
GO  TO  210 
DO  L' 1 . NEON 

Y  f  L  )  -  P  (  L)+TEMP1*<  YP(  I.  )  PH  I  I  L  ,  1  '  > 

END  DO 

CALL  FIX.Y.YP) 

DO  L- 1 , NEQN 

PHI(L,KPl)-YPfL)  PHI(L,  l) 

P  H I (  L,KP2)=PHII  L,KP1)-PHI( L.KP2  ' 

END  DO 
DO  I  -  l  ,  K 

DO  1.  1  ,  NEON 

PH  I  <  L  ,  n  =  PH  I  L  ,  I  '  -PH  I  <  I. ,  KP  1  ) 

END  DO 
END  DO 
ERKP1  0 . 

IE  KNEW  .EQ.  KMl  .OH.  K  EG.  12  PHASE  1  -  .  FALSE . 

IE  ( PHASE  1;  GO  TO  220 
IE  KNEW  . EQ.  KMl)  GO  TO  220 
IE  <  KP 1  . GT .  NS  )  GO  TO  240 

DO  I.  1  ,  NEQN 

ERKP1  ERKP  1  *  PH  I  L  ,  KP2  '  W  T  <  I.  **2 
END  DO 

ERKP1  * ABSH*GSTR' KP1  )  *SQRT  ERKP!  t 
IE  K  L  E .  1  )  THEN 

IE  ERKP  l  .  GE.  .  5  *  E  R  K  1  GO  TO  240 
GO  TO  220 
END  TE 

IE  i  E  R  K  M 1  .  LE  M I N  E R K , E R K P 1  ’  GO  TO  220 
IE  ERKP1  .GE.  EHK  .OR.  K  . EQ .  12  GO  TO  240 

K  KP1 
ERK  ERKP I 
GO  TO  240 
K  KM  1 
ERK - ERKM 1 
HNEW  H  *  H 

IE  PHASE  1)  GO  TO  250 

IE  (P5EPS  GE.  E  R  K  *  T  W  0  K* I  GO  TO  250 

HNEW  H 

IE  < P  5  E  P  S  GE.  ERK  GO  TO  250 

TEMPI’  K*  1 

R-  f  P5EPS / ERK  «  *  1-  TEMP 2 

HN  EW - A8S  H  *M A  X  .500, MINI. 9D0.R 
HNEW-SIGN.  MAX  HNEW,  fOURUAABS  \  ,H> 

H  =  HNEW 
RETURN 
END 

SUBROUTINE  INTRP  X  .  Y  ,  XOIiT  .  YOUT  ,  Y  POUT  ,  NEON  ,  KO  LI)  ,  PH  I  ,  PS  I 
SUBROUTINE  STEP  APPROXIMATES  THE  SOLUTION  NEAR  X  BY  A 
POLYNOMIAL.  SUBROUTINE  INTRP  APPROXIMATES  THE  SOLUTION 
AT  XOUT  BY  EVALUATNG  THE  POLYNOMIAL  THERE. 

IMPLICIT  DOUBLE  PRECISION  A  H,0  7. 

DIMENSION  Y (  NEQN  , YOUT I NEQN  . YPOUTi  NEQN  >  , PHI  NEQN,  10  i  , 


fcPS I ( 12) ,  G  ( 13) ,Wf 13) ,  RHO  ( 13) 

KM A  PHI, XOUT , YPOUT 
G  (  1 )  =  1  • 

RHO ( 1  )  =  1  . 

H I  =  XOUT -  X 
K I =  KOLD  +  1 
KIPl=KI+l 
DO  1=1, KI 
TEMP  1  =  1 
W( I ) = 1 . /TEMPI 
END  DO 
TERM=0 . 

DO  J=2,KI 
JM  1  =  J  -  1 

PSIJM1-PSI(JM1  ) 

GAMMA- ( HI  +  TERM) / PS  I JM 1 
ET  A  =  H I / PS  I JM 1 
LIMIT1-KIP1-J 
DO  I  -  1 , LIMIT1 

W  (  I  ) -G AMMA*W  <  I  1-ETA«W(  1  +  1  ) 

END  DO 

G( J ) =W( 1  ) 

RHO l J  )  -GAMMA *RHO( JMi  ) 

TERM: PS  I JM1 
END  DO 
DO  l- 1 . NEON 

YPOUT ( L  l  : 0 . 

YOUT ( L )  -  0 . 

END  DO 
DO  J  =  1 , K  I 
I-KIP1-J 
TEMP2  T  G (  I  ) 

TEMP3  -  RHO  <  I  ) 

DO  L -  1 , NEON 

YOUT  r  I.  i  -  YOUTf  L  i  ♦TEMP2*PHI  (  L,  I) 

YPOUT  1.  )  =  YPOUT  (  L  )  ♦TEMP3*PHI(  L,  I  ) 
END  DO 
END  DO 
DO  L -  1 , NEON 

YOUT  (  L  )  -  Y  f  l )  *H  I  *  YOUT  (  I,  > 

END  DO 
RETURN 


APPKNDIX  B 


TWriD  I'SKK'S  WAM  AL 


1  .  !nt  rivluct  ton 

rv-l  1  [)  Is  a  one-d  i  mens  i  ona  1  larm-s  Igna  !  1  vs  ;  s  pruTim  f  >  r  helix 


travel! ng-wa'. 

•e  t  uhes  based. 

on  the  theorv 

Vve  1  iped  *■  v  1 

V 

r  We  i 

:id  V. 

K.  Sc  her ha. 

The  p  r  >g  r  im 

can  be  us ed  to 

pred.  let  the  | 

r  f 

t  . i  i:n> 

ot  a 

helix  travel 

I  u  g  —  w  a  v  e  tube 

before  it  is 

is  se’it  ;  ,  d  .  T 

1  !« 

i-'l'.t  t 

■  the 

pr  >gran  cons  i  sts  it  i  n  t .  <  r  :tu-»  t  i  •  -n  ah-mt  the  bear...  th<  i  -.pi.t  s  i  Is  ,  and 
the  helix.  hinnett-rs  Jes,  t  (Mi«  the  helix  ate  t  j  •  i .  . « I  1  -  -Ft  ai  rtvl  1  r  .ra 
1*0 1  d  c  i  r  cu  i  t  testing,  while  intorn.tr  l.in  ah  the  b*  am  -.iv  c  r,e  !i  .m 

he  in  analyzer  results  or  computer  design  data.  i  ri-a  to  running  TV!  in, 
the  engineer  creates  an  input  data  tile  which  lists  these  par  met  ers . 
This  manual  describes  how  tills  tile  is  to  he  created  and  hew  to  run 
TWT1D.  Complete  document  at  i  on  for  t  tie  program.  Including  a  s  ur.e  cade 
listing,  is  given  in  A  FT  HR  report  No.  1',  available  from  t  tie  Flectrieal 
'••Engineering  Department  at  the  Pnlversltv  ot  it  ah. 

2 .  Input  Data  File 

The  input  file  should  he  organized  along  tin-  lines  of  fig.  A.l. 
The  number  of  lines  in  the  file  is  variable  and  is  determined  bv  the 
number  of  frequencies,  velocity  steps,  attenuators,  etc.,  present  in  the 
tube  being  modeled.  The  kev  following  Fig.  A.l  describes  eat  h  if  the 
variables  listed  In  the  input  file,  and  Fig.  A. 2  shows  a  sample  input 
file.  Since  the  input  flies  can  become  lengthy,  it  is  strongly  recom¬ 
mended  that  before  using  TV*T1D,  the  input  data  file  be  checked  with  the 


program  TWTIS.  TWT I N  is_  a  short  program  which  roads  TUT1I) 
and  lists  the  parameters  so  that  the  engineer  can  verity  the 
used.  The  sample  input  tile  shown  In  Fig.  A. 3  was  checked 
and  the  result  Is  shown  in  Fig.  A. a. 


,  r\e 

D  a  t  n 

1 

VO  10 

RAD  Fill  TP  I  ZF 

o 

STMF 

NEMO  N'i  VIA  NUMR 

3 

F  1  1 

p  X  N  1  l  ■  P  H  V  E  L  (  1  ■  A  X  I  M  P  ■  1  R  F  L  S  S  i  1 

4 

ro 

P  I  N  12 

5 

FIT' 

P  I  N  I  .  3 

h 

z  c  l  ; 

C  K  T  R  A  D  1  )  TAPER;! 

7 

ZC  2 

CKTRAP  >  2 )  TAPER ( 2  - 

1  1 

ZC ■ 6 )  CKTHAD i 6 ) 

TAPE  R  *  6  i 

3  ) . 

12 

P  H  V  E  L  N  i  1  ,  1  P  H  V  E  L  N  (  1  ,  2  )  P  H  V  E  L  N  <  1  , 

1  3 

AXIMPNtl.l  AXIMPNt 1,2)  AXIMPNil, 

3  ! . 

1  4 

RFLSSNtl.l  RFLSSN (  1 , 2  >  RFLSSNtl, 

3  ) . 

15 

Z  A  T  N  S  (  1  ■  Z  A  T  N  E  i 

1  )  LAM A  X i  1  ) 

A  T  Y  P  E  ( 

1  > 

1  6 

Z  A  T  N  S  !  2  Z  A  T  N  E  < 

2  )  LAMA  X (  2  ) 

ATYPE ( 

2  ) 

20 

Z  A  T  N  S  6  ■  Z  A  T  N  E  t 

6 )  LAMAX i 6  ) 

ATYPEi 

6  ) 

2  i 

Z  R  S  T  1  Z  R  S  T  (  2  ) 

Z  R  S  T  (  3  )  .  .  . 

. . ZRST ( 

)  ') 

SCS  N E  LEFT  I) PC 

NPR  INTSKP 

RSTART 

FDUMMY 

input  files 
i  n pu  t  he  i  rig 
with  TWT I N , 


THETA-  1 


p  h  v  e  i.  n  (  6 ,  :n 

AX IMPN i 6 , 3 ! 
R  F  L  S  S  N  <  6,3) 

10  ) 


Fig.  A . 1 .  Input  data  file  organization. 


KEY  TO  INPUT  DATA 


Variable 

VO 

10 

RAD 

FILL 

TPI 

ZF 

NUMF 

NUMC 

NUMA 

NUMR 

FI(  i ) 

PINKi) 

PHVEL(i) 

AXIMP(i) 

RFLSS(l) 

THF.TA(i) 

ZC(1) 

CKTRAD(i) 

TAPER!  I  ) 

DHVELN(  1  ,  j  ) 


Desert  pt Ion 


Dc  beam  voltage  (volts) 

Dc  beam  current  (amps) 

Mean  helix  radius  (in) 

Beam  fill  factor 
Helix  pitch  (turns/in) 

Length  of  circuit  (in) 

Number  of  input  frequencies  (3  maximum) 

Number  of  circuit  changes  and  velocity  tapers 
(6  maximum) 

Number  of  attenuators  (6  maximum) 

Number  of  restart  data  printouts  (10  maxi¬ 
mum).  The  restart  printout  contains  all  the 
information  that  is  needed  to  restart  a  run. 

Frequency  of  ith  signal  (GHz) 

Input  power  of  ith  signal  (dBm) 

Normalized  phase  velocity  of  1th  signal  (v  /c) 

Axial  coupling  Impedance  of  ith  signal  (ohms) 

RF  circuit  loss  for  ith  signal  (dB/in) 

Injection  angle  of  ith  signal  (degrees) 

Pos ition  of  ith  velocity  step  or  taper  (in) 

Radius  of  circuit  beyond  ith  change  (in) 

Length  of  velocity  taper  (in).  Use  0  for 

velocity  step. 

Normalized  phase  velocity  of  TPI  of  the  helix 
for  jth  signal  beyond  ith  change.  If  the  TPI 
is  used,  the  characteristics  of  the  new  sec¬ 
tion  are  calculated  by  scaling  the  phase 


velocity  and  impedance  values  from  the  previ- 


ous  section.  When  the  TPI  is  specified, 

use 

J 

-  -j 

an  axial  impedance  of  99.  TPI  scaling 

can 

* . . 
.-.-j 

only  be  used  when  the  TPI  of  the  previous 

section  is  given.  Also,  TPI  scaling  will 

not 

work  when  going  from  a  vaned  to  an  unvaned 

•  1 

circuit,  or  vice  versa. 

•'  '"I 

fjy 

13 

AXIMPN(i.j) 

Axial  Impedance  for  jth  signal  beyond 

ith 

change  (ohms) 

i 

14 

RFLSSN(i.j) 

RF  circuit  loss  for  jth  signal  beyond 

ith 

•*, « 

change  (dB/in) 

;-.c- 

,.v 

15-20 

ZATNS(i) 

Start  position  of  ith  attenuator  (in) 

.‘.■c 

%■ 

ZATNE(i) 

Stop  position  of  ith  attenuator  (in) 

LAMAX(i) 

Maximum  loss  of  ith  attenuator  (dB/in) 

ATYPE(i) 

Type  of  ith  attenuator  (1,  2,  3,  or  4). 

See 

■"v; 

Fig.  A. 2  for  the  different  attenuator  shapes 

available. 

21 

ZRST(i) 

Position  of  1th  restart  data  printout  (in) 

’  ! 

■4 

22 

SCS 

Space-charge  switch: 

• 

0  =  no  space  charge 

1  =  full  space  charge 

2  =  polynomial  approximation  to  space 

charge 


The  space-charge  switch  determines  how  the 
repulsive  force  between  electrons  is  calcu¬ 
lated.  If  the  tube  has  low  beam  current 
density,  or  quantitative  results  are  not 
required,  then  it  is  appropriate  to  ignore  the 
space  charge  and  set  SCS  to  0.  If  the  tube 
has  few  circuit  changes  and  a  small  step  size, 
the  polynomial  approximation  to  the  space 
charge  will  give  reasonable  results,  and  SCS 
should  be  set  to  2.  For  highest  accuracy,  SCS 
should  be  set  to  1  so  Rowe’s  full  space-charge 
expression  is  used. 


NELECT  Number  of  electron  disks  (16,  32,  or  64).  Use 

16  for  one  input  signal;  '32  for  two  input 
signals,  and  64  for  three  input  signals. 


DPC 


Printout  interval  (in) 


% 


-  73  - 


NPR 

INTSKP 


RSTART 


FDUMMY 


Printout  option  1,  2,  3,  or  4  (see  Section  IV 
for  description) 

Number  of  skips  between  evaluation  of  space- 
charge  function.  Since  the  electron  positions 
do  not  vary  rapidly,  the  space-charge  function 
can  be  evaluated  intermittently  to  reduce 
computation  time. 

Restart  data  input  switch: 

0  =  no  restart  data 
1  =  read  restart  data 

Set  to  1  if  the  run  is  to  be  restarted.  A 
second  data  set  must  be  included  after  the 
first  one. 

Dummy  frequency  (GHz).  Any  data  set,  which 
specifies  multiple  frequencies  which  are  not 
harmonically  related,  must  also  specify  a 
dummy  frequency.  This  dummy  frequency  is  the 
highest  frequency  of  which  all  Input  signals 
are  harmonics.  The  lower  this  frequency,  the 
less  accurate  the  results  will  be.  The  dummy 
frequency  is  not  followed,  as  it  is  an  artifi¬ 
cial  construction. 
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ATYPE  3  Output  attenuator  ATYPE  4  Uniform  attenuator 


Fig.  A. 2.  Attenuator  shapes 


10200.  .280  .0379 

.5  19.18  9.5 

13  11 

18.  10 .  . 1864  3 . 6 

.75  0  . 

4.68  .0379  0 

5.30  .0379  0 

8.00  .0379  .5 

3.75  3.44  3.00 

3.00  3.50  100  2 

4 . 5 

116.212000 

Fig.  A. 3.  Sample  input  file 


INPUT  DATA  FOR  FILE  SAMPLE.  IN 


HELIX  VOLTAGE! Volts)- 
CATHODE  CURRENT ( Amps ) = 
MEAN  HELIX  RADIUS! in) - 
BEAM  FILL  FACTOR  : 
HELIX  PITCH(tpi)= 
LENGTH  OF  CIRCUIT! in ) - 
NUMBER  OF  FREQUENCIES^ 
NUMBER  OF  CHANGES = 
NUMBER  OF  ATTENUATORS  = 
NUMBER  OF  RESTART  DATA  PRINTOUTS^ 
NUMBER  OF  ELECTRON  DISKS= 
SPACE  CHARGE  SWITCH- 
CIRCUIT  PRINT  INTERVAL!  in = 
PRINT  OUT  OPT  I  ON  = 
NUMBER  OF  SPACE  CHARGE  SKIPS= 


10200 . 00 
.  2800 
. 03790 
.  5000 

i  9 .  180 

9 . 5000 
1 


FREQ. 
!  GHz  > 
18.00 

CHANGE 


INPUT  PWR . 
(  dBm  ) 

10 . 000 


PHASE  VEL. 

<  Vp/c ) 

.  1864 


AXIAL  IMP. 
( ohms ) 
3.600 


CKT.  LOSS 
( dB/  in ) 

.  750 


REI..  PHASE 
( degrees ) 

0 . 000 


POSITION  CKT.  RADIUS  TAPER  LENGTH 


CHANGE 


ATTEN  . 


RESTART 


(in) 

(in) 

>in: 

4 . 680 

.0379 

0  .  0 

5 . 300 

.  0379 

0  .  0 

8.000 

.0379 

.  5 

PHASE  VEL. 

AXIAL  IMP. 

CKT.  LOSS 

( Vp/c  ) 

( ohms  ) 

(  d  B  /  l  n  ' 

.  1988 

3 . 75 

.  7  5 

.  1796 

3  44 

.  7  5 

.  1827 

3  .  00 

.  7  5 

:art  stop 

MAX.  LOSS 

TYPE 

in)  (in) 

(dH  ) 

000  3.50 

100.0 

n 

POSITION 

(in) 

4 . 5000 


Fig.  A. 4.  TWTIN  run  for  sample  input  file 


3.  Running  the  Program 


Once  an  Input  file  has  been  created  and  checked,  the  engineer  is 


ready  to  run  TWT1D.  When  run,  the  program  prompts  the  user  for  the  name 


of  the  input  data  file,  the  output  data  file,  and  the  plot  file.  The 


output  data  file  is  where  all  the  data  from  the  run  will  be  stored,  and 


the  plot  file  will  contain  the  data  needed  to  make  plots  of  power  versus 


position,  and  phase  angle  versus  position.  Although  any  name  can  be 


used  for  the  output  and  plot  files,  the  user  is  cautioned  against  using 


names  of  existing  files.  That  can  result  in  the  original  contents  of 


those  files  being  lost.  It  is  recommended  that  the  names  of  the  output 


and  plot  files  be  related  to  each  other  so  that  it  is  easy  to  recall 


which  plot  file  corresponds  to  which  output  file.  For  example,  if  the 


input  file  is  called  RUN1.DAT,  the  output  and  plot  files  could  be  named 


RUNl.OUT  and  RL'Nl.PLT,  respectively.  Figure  A. 5  shows  the  terminal 


display  for  a  TWT1D  run.  As  the  program  runs,  the  axial  position  and 


power  will  be  printed  to  the  terminal. 


TWT1D  10:50  AM  MON.,  20  JAN,  1986 

INPUT  DATA  FILE  NAME?  SAMPLE. IN 

OUTPUT  FILE  NAME?  SAMPLE. OUT 

POWER  PLOT  FILE  NAME?  SAMPLE. PLT 

ALL  OUTPUT  DATA  WILL  BE  ROUTED  TO  FILE  SAMPLE. OUT 
ALL  PLOT  DATA  WILL  BE  ROUTED  TO  FILE  SAMPLE. PLT 


Fig.  A. 5.  Terminal  display  for  a  TWT1D  run. 

(User's  input  is  underlined.) 

4 •  Output  Data  File 

The  TWT1D  output  file  for  the  sample  input  file  is  shown  in  Fig. 
A. 6.  Notice  that  at  each  print  interval,  several  variables  describing 
the  performance  of  the  TWT  are  printed  to  the  output  file.  Which  vari¬ 
ables  are  printed  is  determined  by  the  variable,  NPR,  in  the  input  data 
file.  The  four  options  available  are  described  and  illustrated  in  this 


section. 


TWT1D  10:58  AM  MON.,  20  JAN.,  1'JHb 
RESULTS  FOR  INPUT  FILE  SAMPLE.  IN 


FREQ. 

18. 


PLASMA 


*Hz)  INPUT  PWR . !  dBm )  REL 

1000  10.0000 

CATHODE  CURRENT! Amps) 
HELIX  VOLTAGE! Volts) 
EFFECTIVE  VO LT AGE ( V o 1 t s  ) 
NORMALIZED  BEAM  VELOCITY 
HELIX  PITCH!tpi 
PLASMA  FREQUENCY (raas/s) 
FREQUENCY  REDUCTION  FACTOR 
SPACE  CHARGE  PARAMETER 
MEAN  HELIX  RADIUS(in) 
BEAM  FILL  FACTOR 
CIRCUIT  PRINT  INTERVAL! in ) 
NUMBER  OF  ELECTRON  DISKS 
NUMBER  OF  ATTENUATORS 
NUMBER  OF  CHANGES 


PHASE! rad ) 

0 . 0000 

.  2800 
10200 .00 
9808.02 
.  1959 
19.180 

.  1 141E  +  1  1 
.  5987 
.  7993 
. 03790 
.  5000 
.2000 


N  NORM. 
1 


PHASE  VEL. (Vp/c) 
.  1864 

PDB  B 

9.874  1.60 

9.909  1.60 

_J0. 32>v 


RESTART  INFORMATION  FOR  Y- 


AX  I AL  IMP . ( ohms ) 
3.6000 


GAIN  PARAMETER 
.0318 


THETA 

-.505 

-1.035 


7.009 


DPDB 
-.628 
.  176 


DTDZ 

-1.0610 

-1.1227 

-kJ4i2- 


ETA 
.  0000 
.  0000 
.  0  00  0- 


. 7872E-0 1 


DA  THETA 

.  3605E-0  1  - . 1 192E  +  02 


DTHETA 
1 975  E  0  1 


ELECTRON  PHASE  POSITIONS 
.  1200E  +  02 
.  1241E  +  02 
.  1280E+02 
.  1318E  +  02 
. 1 353E+02 
.  1  388E  +  02 
.  1422E  +  02 
.  1458E+02 
.  1494E+02 
.  1532E  +  02 
.  1572E+02 
. 1614E+02 
.  1656E  +  02 
.  1700E+02 
.  1743E  +  02 
.  1 7  86E  +  02 

ELECTRON  VELOCITIES  „„ 

. 7256E-0  1  . 4444E-0  I  .6205E-02 

-.7122E-01  -.9445E-01  -.1012E+00 

-.  7145E-01  - . 4 34 3E - 0  1  -.1196E-01 

.4919E-0  1  .  7247E-0  1  .8605E-01 


4 . 680 
4.880 


9.350 

9.500 


32-075 

33.921 


542127 

54.201 


1.60  -12.394 

-  .45  -12.335 


2.41  -28.754 

2.41  -29.300 


3519E-01 
9222E-0 1 
1 990E-0  1 
.  8667E-0 1 

9.327 
9. 232 


1.430 
.  497 


•1  .  1459 
-  .  4336 


1  . 9667 
■  1  . 2292 


.  0006 
.0009 


>tr906 

.0921 


ENDING  PRINT  OUT 
SATURATED  POWER  (dBm)= 
GAIN  AT  SATURATION  (dB>  = 
SATURATION  EFFICIENCY  (*)= 
LENGTH  AT  SATURATION  (inl= 
SMALL  SIGNAL  GAIN  !dB)  = 

TWT1D  1:35  PM  MON. ,  20 


Fig.  A. 6.  TV/T1D  output  file. 
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This  Is  the  printout  option  which  will  most  often  be  used  by  tube 


designers.  It  prints  out  the  following  variables. 

Z  Axial  position  (in) 

PDB  Power  (dBm) 

B  Pierce's  velocity  parameter 

THETA  The  phase  angle  of  the  circuit  wave  (radians) 

DPDB  The  partial  derivative  of  power  with  respect  to  position 

DTDZ  The  partial  derivative  of  phase  angle  with  respect  to 
position 

ETA  The  efficiency 

D  Pierce's  loss  parameter 

NPR  -  2 

This  option  prints  out  all  the  variables  in  NPR  =  1 ,  as  well  as 
information  about  energy  in  the  beam  and  in  the  circuit  wave.  This  can 
help  indicate  whether  or  not  the  computer  simulation  is  a  good  repre¬ 
sentation  of  the  physical  process.  The  variables  are: 

EBEAM  Power  lost  by  the  beam,  normalized  by  the  total  beam 

power 

EWAVE  Power  contained  in  the  circuit  wave  and  in  the  field 

between  the  circuit  and  the  beam 

EBLNC  (IeBEAm|  -  | EWAVE | ) / ( | E BEAM |  +  EWAVE|) 

If  the  simulation  is  reasonable,  EBLNC  will  be  near 
zero.  A  "good"  value  for  EBLNC  is  on  the  order  of  0.1. 
However,  this  parameter  does  not  account  for  the  circuit 
wave  power  lost  in  attenuators  so,  in  those  regions, 
EBLNC  will  not  be  near  zero. 


a 


NPR  -  3 


This  option  Includes  all  the  variables  in  NPR  •  1,  as  well  as 
Information  for  creating  electron  flight  diagrams  and  Cutler  charts. 
The  variables  are: 

PHI(t,L)  Electron  phases  for  the  fundamental  frequency 

2*C(l)*U(i)  Normalized  electron  velocities 

NPR  -  4 

This  option  will  print  out  all  the  variables  listed  above.  It 


generates  a  lot  of  data!  All  of  the  NPR  options  are  illustrated  in  Fig. 


/ 

PDB 

R 

THE  T  A 

.ip I)  H 

i)  T  D  / 

ETA 

.  200 

9.874 

1  .  60 

6  0  6 

.  6  2  8 

1  .0610 

.  0  0  0  0 

4  00 

9 . 90  9 

1  .  60 

1  .  U  <6 

.  1  76 

1  .  1227 

uOOO 

PR  2 

7 

PDB 

B 

THETA 

DPI'R 

DTI)/ 

ETA 

ERE  AM 

EWA  VE 

KB  I  NC 

.  200 

9 . 874 

1  .  60 

.505 

.628 

1  .0610 

.0000 

■  .  1  0  2  4  5  E 

07  -  .  10897E  06 

-  8281 1 

.  400 

9.909 

1  .  60 

1.045 

.  i  7  8 

1  .  1227 

.  0000 

.  14947E  06  -  .  90  00 3 E  07  .  24522 


7 

PDB 

B  THETA 

DRUB 

DTDZ  ETA 

PHI 

I , 1  ELECTRON 

PHASE  ARRAY 

t  FUNDAMENT  AL  1 

2  *C 

l  *  i :  I  ; 

NORM  . 

AC  VELOCITY 

ARRAY 

.200  9 

.874 

1  .  60 

.  5  0  5 

.628 

1.0610  .0000 

.  5049 

.  8973 

1  .  290 

1  .  682 

2.075 

2.468 

2.861 

3  .  254 

9  .  64  7 

4.040 

4.44  3 

4 . 825 

5.218 

5.611 

6.003 

6. 396 

. 4965E  04 

.  9  1  3  3  E  -  0  4 

.  1291 E  0  3 

-  .  1  4  7  2  E  -  0  3 

.  1  4  30  E  03 

- . 1 169E-03 

.  7307E- 04 

- . 1809E-04 

. 3  9  6  5  E  04 

. 9 1 3  4  E  0  4 

.  1291E-03 

.  1472E-03 

. 1429E-03 

. 1 169E-03 

.  7  3  0  2  E  -  0  4 

.  1  806E-04 

PR 

4 

Z 

PDB 

B  THETA 

D  PD  B 

DTDZ  ETA 

EBE  AM 

EWAVE  EBLNC 

PH  I 

(1,1-  ELECTRON 

PHASE  ARRAY 

( FUNDAMENTAL) 

2  *  C  f  1  '  *  U  f  I  ) 

NORM. 

AC  VELOCITY 

ARRAY 

.  200  9 

.  874 

1.60 

.  505 

-  .  628 

-1.0610  .0000 

10245E-07 

10897E-06 

-.  8281  1 

.  5049 

.  8973 

1  .  290 

1 . 682 

2.075 

2 . 468 

2.861 

3.254 

3.647 

4 . 040 

4.433 

4.825 

5.218 

5.611 

6.003 

6.396 

- 

. 3965E-04 

-. 9133E-04 

-  .  1291E-03 

1472E-03 

- 

.  1430E-03 

- . 1 169E-03 

- . 7307E-04 

1809E-04 

. 3965E  04 

. 9134E-04 

. 1291E-03 

.  1472E-03 

.  1429E-03 

. 1 169E-03 

. 7302E-04 

.  1806E-04 

Res'  » r  t  i  ng  a  Run 

i  r  ten  an  ei^l:;eer  will  want  t examine  the  effect  of  different 
•  tput  helices  an  a  tube's  performance.  TWT1D  allows  the  user  to  do 
lets  wit". ait  having  t  >  rerun  the  Input  helix  In  each  case. 

T*  •  a  run  Is  to  he  restarted,  the  input  file  for  the  initial  run 
s’.oul  1  have  the  number  of  restarts,  Nl'MR  set  to  1,  and  the  position  of 
'  he  restart  printout,  ZRST,  set  to  the  position  of  the  restart.  The 

out  put  file  f  or  that  run  will  contain  all  the  Information  needed  to 

restart  the  run.  A  restart  printout  at  4.5  inches  is  illustrated  in 

Fig*  A.b.  To  restart  the  run,  the  data  from  the  restart  must  be 
appended  to  the  original  input  file,  and  the  restart  data  input  switch, 
RSTART,  set  to  1.  Figure  A. 8  shows  the  input  file  for  a  restarted 
run.  Note  that  all  the  character  data  have  been  removed  from  the 


restart  data. 
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